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Abstract. This paper is devoted to the asymptotic analysis of boundary layers in 
periodic homogenization. We investigate the behaviour of the boundary layer corrector, 
defined in the half-space fl n , a '■= {y • n — a > 0}, far away from the boundary and prove 
the convergence towards a constant vector field, the boundary layer tail. This problem 
happens to depend strongly on the way the boundary dQ. n , a intersects the underlying 
microstructure. Our study complements the previous results obtained on the one hand 
for n £ RQ d , and on the other hand for n ^ EQ d satisfying a small divisors assumption. 
We tackle the case of arbitrary n £ RQ d using ergodicity of the boundary layer along 
dQ. n ,a- Moreover, we get an asymptotic expansion of Poisson's kernel P — P(y,y), 
associated to the elliptic operator —V • A(y)V- and 0„ >a , for \y — y\ — > oo. Finally, we 
show that, in general, convergence towards the boundary layer tail can be arbitrarily 
slow, which makes the general case very different from the rational or the small divisors 
one. 



In this paper we investigate the behaviour of vu = vu{y) G R , solving the elliptic 
system 



with periodically oscillating coefficients and Dirichlet data, far away from the boundary 
of the half-space {y ■ n — a > 0} C R d . Understanding these asymptotics is an important 
issue in the study of boundary layer correctors in periodic homogenization. 

Throughout this paper, N > 1, d > 2, n € S d_1 is a given unit vector and a G R. 
The Dirichlet data vq = Vo(y) G M. N is defined for y G M. d ; so is the family of matrices 
A = A a P(y) G Mn (R) indexed by 1 < a, (3 < d. Small greek letters like a, (3, 7, rj usually 
denote integers belonging to {1, . . . d}, whereas i, j, k stand for integers in {1, . . . N}. 
Therefore, taking advantage of Einstein's convention: 



The main assumptions on A are now 

(Al) ellipticity: there exists A > such that for every family £ = £ a G R N indexed 
by 1 < a < d, for all y G R d , 



1. Introduction 



(1.1) 




v b i =0, y • n - a > 
vu =vo(y), y-n-a = 



V ■ A{y)Vv h ^= d ya (Af{y)d yp v bhj 



(A2) periodicity: A is 1-periodic i.e. for all y G R d , for all £ G Z d , 



A(y + = A(y); 



(A3) regularity: A is supposed to belong to C°° (R d ) . 



Moreover, we assume 

(Bl) periodicity: vq is a 1-periodic function; 
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(B2) regularity: vq is smooth. 



Before focusing on system (1.1) itself, let us recall the main steps of the homogenization 
procedure leading to the study of Vu- We carry out our analysis on the elliptic system with 
Dirichlet boundary conditions 

-V • A (§) Vu e = f, xen 

u e = o, x £ on 



(1.2) 



posed in a bounded Lipschitz domain Q, C M. d endowed with a periodic microstructure. 
This system may model, for instance, heat conduction (JV = 1), or linear elasticity (d = 2 
or 3 and N = d) in a periodic composite medium. We are interested in the asymptotical 
behaviour of u £ , when e — > 0. This problem is of particular importance when it comes 



to the numerical analysis of (1.2). The oscillations at scale e are too fast to be captured 
by classical methods. One therefore looks for an approximation of u £ , taking into account 
the small scales, without solving them explicitely. To gain an insight into these numerical 
issues, we refer to |SV06| ; for a general overview of the homogenization theory see |BLP78j 
or |CD99| . 

1.1. Some error estimates in homogenization. For a given source term / G L 2 (S1), 
boundedness of and ellipticity of A yield the existence and uniqueness of a weak solution 



u 6 in Hq(Q) to (1.2). Moreover, we deduce from the a priori bound 

I« £ IIh-i^ < c 



l-HS(fi) - ° WJ IIl 2 (H) i 

where C > is a constant independent from e, that up to the extraction of a subsequence, 
u 6 converges weakly in Hq(Q). A classical method to investigate this convergence is to 
have recourse to multiscale asymptotic expansions. We expand, at least formally, u £ in 
powers of e 

(1.3) u £ rj u° (x, J) + eu 1 (x, ^) + e 2 u 2 (x,^j+... 

assuming that u l = u l (x, y) G ~R N is 1-periodic with respect to the fast variable y. Plugging 



(1.3) into (1.2 ), identifying the powers of e and solving the cascade of equations then gives: 
(1) that u° does not depend on y and that it solves the homogenized system 

-V • A°X7u° = /, ieS! 

u° =0, x£dn 

where the constant homogenized tensor A = A°' a ^ G (R) is given by 
A °,*P := [ A^(y)dy + [ A^(y)d y y(y)dy, 



and x = X^iu) £ Mjy (^) is the family, indexed by j3 G {1, . . . d}, of solutions to 
the cell problem 

rd 



n4) { -V y -A(y)V yX P =d ya A^ , yGT d 

1 ' I W(v)«fo =0 



(2) that for all x G Q and y G T d , ^{x^y) = X a (y)d Xa u°(x) + u l (x); 

(3) and that for all x G 0, and y G T d , u 2 {x,y) := T a ^ {y)d Xa d xp u {x)+ X a {y)d Xa u 1 {x) + 
u 2 (x), where T = T a/3 (y) G M N (M) is the family, indexed by a, /3 G {1, . . . (i}, 
solving 

-VyA(y)V y T^ =B^-J Jd B^(y)dy, y G T d 
f Td r a P(y)dy =0 

and 

B°^(y) := A^(y) + ^(yH^y) + ^ (^ a (y)/(y)) . 
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-2,e 



u £ (x)—u°(x)—ex (§) Vit°(a;) 



One can then carry out energy estimates on the error r bl 
£ 2 T (|) • V 2 it°(x) solving the elliptic system 

-V • A (§) Vr 2 < £ = f e , xen 

r 2,e _ g£^ x E 8U ' 

To bound r 2,e one needs some regularity on u°, say u° £ H 4 (£l) (for refined estimates, 
involving lower regularity on u°, see [MV97| iPrallj ). Under this coarse assumption 



11/ 



E V ( n) < Ce \\u°\\ H i {n) 



and 



\Q i 



X 



■VV 



< C£2 



.01 



One can therefore show that 
(1.5) 



£X 



< 



< C£2 



H2(n) 11 H# 4 W 

, which implies 



Jf!(f2) 





■ vV 







H!(n) 



H 4 (n) 



and 



lL 2 (n) 



0(ei). 



The latter estimate shows, that the zeroth-order term u° is a correct approximation of u e . 
Estimate (1.5 1 is however limited by the trace on of u 1 (•, -). Actually, the periodicity 
assumption of the ansatz ( 1.3 ) with respect to the microscopic variable y, is not compatible 
with the homogeneous Dirichlet condition on the boundary. In order to force the Dirichlet 



condition one introduces a boundary layer term u bl 
expansion (1.3). More precisely, ulf 

( —X7 . A ( 

(1.6) 



1,5 



1.E 



(x) G 



at order e 1 in the 



solves 
V-A(f)V«£ 



0, 



-u [x 



.1) 



l£0 



and the corrected Ansatz is 



(1.7) 



u 



u° (x) + e 



u x 



+ u 



1,£ 



+ 



Adding the boundary layer at first-order improves (1.5): if u € H (£1), 



U 



ur 



u 



n 



- ] • Vn° 



eu 



l,e 
6/ 



However, such a trick remains useless as long as one is not able to describe the asymptotics 
of u^f when e — > oo. 



1.2. Homogenization of boundary layer systems. As for (1.2), the problem is to show 
that (1.6) can be in some sense homogenized. A few remarks are in order: 

(1) System (1.6) exhibits oscillations in the coefficients as well as on the boundary. 

(2) The oscillations along <9f2 are not periodic in general. This can be expressed by 
the fact that the boundary breaks the periodic microstructure. 



(3) The a priori bound on u,l 



u 



l,e 



fl"!(n) 



< C 



U X 



0(e-h) 



does not provide a uniform bound in H 1 ^). 



These issues make the homogenization of the boundary layer system ( 1.6 ) far more difficult 



than the homogenization of (1.2). The results obtained in this direction are still partial. 

.rst st( 
in e. If N 



The first step towards the asymptotic behaviour of ulf is to get a priori bounds uniform 



1, the maximum principle furnishes a uniform bound in L°°(Q) on u^f. 



A way of investigating (1.6) in the case when N > 1 is to represent uh £ in terms of 



the oscillating Poisson kernel associated to —V • A (|) V- and Q. In the series of papers 
[AL87, AL89] IAL91] . Avellaneda and Lin manage to get estimates, uniform in s, on these 
Green and Poisson kernels, as well as expansions valid in the limit e — > (for recent 
progress in this direction, see [KLS12J). One of the results of [AL87] (see theorem 3) is the 

< C, for 1 < p < oo, valid under the assumption that 80, is 



uniform bound 



l,e| 



\lp(q) 



at least C 1,a , with < a < 1. Note that an L 2 (Q) bound on u^f yields an H 1 ^) bound 
on the gradient, for oj <s= $7 compactly supported in Q (see [AA99J IGVMll] ). The strong 
oscillations of Vu 6 | E are filtered out in the interior of the domain and concentrate near the 
boundary. Hence, the multiscale expansion is right up to the order 1 in e in the interior: 



(1.9) 



< Ce\ 



u 



However, to improve the asymptotics up to the boundary, one needs another approach. 
Namely, we seek after a 2-scale approximation of the boundary layer corrector 



u 



v x 



Take xq G dO, a point at which there exists a tangent hyperplane directed by n := n(xo) G 
S^" 1 and assume that Q is contained in {x ■ n — xq ■ n > 0}. Then plugging formally v into 
(1.6) gives at order e~ 2 



(1.10) 



-V y ■ A(y)V y v(x ,y) 



n-x y Q 





0, y ■ n 

v(x ,y) = -x{v) ■ V«°(x ), yn-^fa 

The variable xq is nothing more than a parameter in this system. If £1 is convex, the 
boundary layer is approximated in the vicinity of each point xq of the boundary by a 
v x = v(xq,-) solving (1.10). This formal idea has been made rigorous for polygonal 
convex bounded domains SlcK 2 , for which only a finite number of correctors v k has to 
be considered, one for each edge K k . By linearity, v k factors into 

v k {x,y) = vlj a {y)d Xa u°{x) 



where for all a = 1, 



(1.11) 



d-> v hi. a = v bi a (y) e ~^- N sorves 



■bl 

-V • A(y)Vv k b f 

k,a 

v bi 



x a (y), 



y ■ n" 
y ■ n k 



a > 
a = 



with a := — . Dropping the exponents k and a, we end up with (1.1). 



System (1.1) is linear elliptic in divergence form. The main source of difficulties one 
encounters is the lack of boundedness of the domain ^l nA '■= {y ■ n — a > 0}. This 
complicates the existence theory, but even more the study of the asymptotical behaviour. 



Moreover, the analysis of ( 1.1 ) depends, in a nontrivial manner, on the interaction between 



dil n ^ a and the underlying lattice. So far, it has been carried out in two different contexts: 
(RAT): the rational case, i.e. n G RQ d ; 

(DIV): the small divisors case, when there exists C, r > such that for all £ G 



Z d \ {0}, for alH = 1, . . . d - 1, 



-d-r 



(i.i2) h-£i>cier 

where (m, . . . nd-i, n) forms an orthogonal basis of . 



Assumption (1.12) means that the distance from every point, except 0, of the lattice Z , to 



the line {An, A G M}, is in some sense bounded from below. Note that this condition, albeit 
generic, in the sense that it is satisfied for almost every n G S d_1 (for more quantitative 
results, see [GVMJ), is not fulfilled by every vector n ^ lR.Q d . 



Let us explain why the description of the asymptotics of Vbi far away from the boundary 
is a crucial step in the homogenization of (1.6). Roughly speaking, one proves in both 



contexts (RAT) and (DIV), that there exits a smooth vu solving (1.1) and that this 
solution converges very fast, when y ■ n 



oo, towards a constant vector field G 



siV 



the boundary layer tail. For a polygonal domain Q with edges satisfying for instance the 
small divisors assumption, we approximate u^f by u 1 solution of 



(1.13) 



Indeed, 



-V • ^Vn 1 = 0, 



u 



l.e 



U 



L 2 (n) 



< 



Hi, 



4 



,k,oo 



,fc,oo 



•Vn l 



xedQnK k 



L 2 (n) 



+ 



4 



,,k,oo 



L 2 (Q) 



and using the decay of the boundary layer correctors in the interior of f2, one proves: 

Theorem 1.1 ( [Prallj . theorem 3.3). If u° G H 2+U) (Q), with to > 0, i/ien i/iere exists 
k > such that 



L 2 (n) 



O (e K ) . 



This shows that the oscillating Dirichlet data of (1.6) can be homogenized. The limit 



system (1.13) involves the tails of the boundary layer correctors. The corrector u 1 has 



been used to implement numerical shemes in |SV()6 . More recently, the obtention of an 
[TTJ-like theor em for a smooth uniformly convex two-dimensional domain has been achieved 
in |GVM| . Again, it relies on the approximation of the domain by polygonals with edges 



satisfying the small divisors condition, emphasizing the key role of systems (1.1). 



We devote section [2] of this manuscript to a review of the previous results obtained on 



(1.1), with an emphasis on the techniques used in the small divisors case. Let us give an 



insight into these theorems: 



(RAT): Among the rich litterature about this case, we refer to [AA99J (lemma 4.4), 
[MV97J (appendix 6) and the references therein. A precise statement is given in 
theorem |2.1 1 It consists of two parts: 

existence: there exists a variational solution vu G C°° (0 n>o ) to (1.1) (unique 

if appropriate decay of Vt>^ is prescribed) ; 
convergence: there exists a constant vector, called boundary layer tail, v^f de- 
pending on a such that vu{y) — v hi°° anci its derivatives tend to exponentially 
fast when y ■ n — >■ oo. 

(DIV): This case was treated in the recent paper |GVMll] by Gerard- Varet and 
Masmoudi. A precise statement is given in theorem 2.2 It consists again of two 
parts: 



existence: there exists a variational solution v^i G C°° (£l n ,a) to (1.1) (unique 

if appropriate decay of Vt> u is prescribed) ; 
convergence: there exists a boundary layer tail independent from a such 

that Vbi{y) — Vjjf and its derivatives tend to when y ■ n — > oo, faster than any 

negative power of y • n. 

The main difference is that the boundary layer tail depends on a in the rational setting 
and not in the small divisors one, which implies that the homogenization theorem 1.1 is 
true up to the extraction of a subsequence e n in the former. In both cases, one can come 
down to some periodic framework to prove the existence of a variational solution. Fast 
convergence follows from a St-Venant estimate. We come back to these points in detail in 
section [2] 

5 



1.3. Outline of our results and strategy. Our goal is to analyse (1.1 ) in the case when 
n ^ MQ d does not meet the small divisors assumption (1.12). Again, one first wonders 



if the system has a solution. This question can be investigated by methods analogous to 
those of |GVMll] , Indeed, their well-posedness result does not rely on the small divisors 
assumption. The latter hypothesis is however essential to show the convergence towards 
the boundary layer tail in the work of Gerard- Varet and Masmoudi. We therefore have 
recourse to another approach based on an integral representation of the variational solution 



to (1.1) by the mean of Poisson's kernel P = P(y,y) £ Mtv(M) associated to —V • A(y)V- 



and the domain £l nA - Basically, 

v u{y) = \ P(y,y)vo(y)dy, 

for every y £ Q n a . At first glance, if n = the d-th vector of the canonical basis of 
if a = and y = (0,e _1 ), e > 0, then 

vu (o,^) = ^0(0,1),^, 0)^v (y,0)dy' 

= L 1 ^ p G (o ' i ^^ o) ) wo G ( "' o) )^ 

Examining the asymptotics of vw far away from the boundary, requires subsequently to 
understand the behaviour of the oscillating kernel ^hrP (f , § ) when e — > 0, or equivalently 
the asymptotical comportment of P(y,y), when y ■ n — > oo and y £ dQ n Q. This is done 
in section [5] relying on ideas and results of Avellaneda and Lin |AL871 IAL91] . We prove 
an expansion for P associated to the domain fi raj o f° r arbitrary n £ S^" 1 . To put it in 
a nutshell, one demonstrates that there exists k > such that for all y £ Sl nj o, for all 
y £ dn n . , 

C 

\P(y,y)-P exp (y,y)\ < ^ _ 

where P exp = Pexp(y,y) is an explicit kernel, with ergodicity properties tangentially to the 
boundary. The precise statement of this key result is postponed to section [5[- see theorem 
\5.3\ We have an explicit expression for the corrector terms. Although this expansion 
is stated (and proved) for the domain f2 raj o it extends to O nia by a simple translation. 
Studying the tail of vm now boils down to examining the limit when y ■ n — > oo of 



[Pex P (y,y)}v (y)dy + / R(y,y)v (y)dy 

with rest R = R(y, y) £ Mjv(K) satisfying \R(y, y)\ < ^ y _y^-i+ K ■ The rest integral tends to 
0. One takes advantage of the oscillations of vq on the boundary to show the convergence 
of the corrector integrals by the mean of an ergodic theorem. Doing so, we demonstrate 
the following theorem, which is the core of this paper: 

Theorem 1.2. Assume that n MQ d . Then, 



[1) there exists a unique solution £ C°° (O na ) D L°° (O njQ ) of (1.1) satisfying 



( L14a ) ll V ^llL-({ r n-t>0}) ^ °. 

( 1 -1 lb) / li<V'M||£-({ 2/ .n-t=0}) dt < OO, 

(2) and a boundary layer tail £ M. , independent from a, such that for all y £ n n>a 



vbi(y) — > v bl 

locally uniformly in the tangential variable. 



Furthermore, one has an explicit expression for «u (see (6.4)/ 



The use of the ergodic theorem to prove the convergence does not yield any rate. One 
wonders therefore how fast convergence of v^i towards vff is. A partial answer is given by 
the next theorem, whose proof is addressed in section [7] 

Theorem 1.3. Assume that d = 2, N = 1, A = I2 and that n ^ IRQ 2 does not satisfy 



(1.12). Then for every I > 0, for all R > 0, there exists a smooth function v$ and a 
sequence (£m)m G]a, oo[ N such that: 

(1) (tM)j[.f i s strictly increasing and tends to 00; 



(2) the unique solution vu of (1.1), in the variational sense, converges towards as 



y ■ n ->■ 00 and for all M G N, for all y' G <9f2 nj0 n B(0, R), 

\vu(y' + t M n) -Vu\ >t]j. 

This result means that convergence can be as slow as we wish in some sense: for fixed 
n ^ Q~R d which does not satisfy the small divisors assumption, for every power function 
1 1 — y t with I > there exists a 1-periodic v$ = vo(y) G M, such that the solution vm of 

-Avu =0, y ■ n — a > 
vm =vo(y), y-n-a = 

converges slower to its tail than {y ■ n)~ l when y-n — > 00. The main obstruction preventing 
Vf,i from converging faster in general lies indeed in the fact that the distance between a 
given point £ G Z d \{0} and the line {An, A G M} is not bounded from below. This is in big 



contrast with the small divisors case, where (1.12 ) asserts the existence of a lower bound. It 



underlines the strong dependence of the boundary layer on the interaction between dQ r 
and Z d . 



1.4. Organization of the paper. In this paper we address the proofs of theorems 



5.3 



1.3 



(expansion of Poisson's kernel), 1.2 (convergence of the boundary layer corrector) and 
(slow convergence) . Section [2] is devoted to a review of the rational and small divisors 
settings. We insist on the existence proof in the non rational case and underline the role 
of the small divisors assumption in the asymptotic analysis. In section |3j some essential 
properties and estimates on Green and Poisson kernels associated to elliptic operators with 



periodic coefficients are recalled. We prove a uniqueness theorem for (1.1), which makes 
it possible to rely on Poisson's formula to represent the variational solution of ( |1.1[ ) and 
explain to what extent the description of the large scale asymptotics of Poisson's kernel 
P boils down to an homogenization problem. The latter is the focus of sections [4] where 
we study a dual problem, and [5] in which an asymptotic expansion of Poisson's kernel, for 



y ■ n — > 00, is established (see theorem 5.3). This work is the central step in our proof 



of theorem 1.2 The last step is done in section [6J where we prove on the one hand the 
convergence towards the boundary layer tail using theorem |5.3| and on the other hand the 
independence of v^i from a. Section [7] is concerned with the proof of theorem 1.3 



1.5. Notations. The following notations apply for the rest of the paper. The half-space 
{y ■ n — a > 0} is denoted by Q n ,a and in the sequel, for any y G £l n ,a and r > 0, 
D(y, r) := B(y, r) n £l n ,a and V(y, r) := B(y, r) Pi d£l n ^ a - The case when a = is frequently 
used: Q n fl =: VL n . For a function H = H(y,y) depending on y, y G M rf we may use 
the following notation: d\ )OL H := d Va H (resp. d2, a H := dy a H) for all a G {1, . . . d}. 
The vectors e\, . . . represent the canonical basis of IR d . The matrix M G Md (K) is an 
orthogonal matrix such that Me^ = n; N G -M^-i is the matrix of the d — 1 first 
columns of M. All along these lines, C > denotes an arbitrary constant independent 
from e. 
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2. Review of the rational and small divisors settings 



In this section we make a short review of the mathematical results on system (1.1). We 



concentrate on giving precise statements for the theorems announced in the introduction 



(cf. section 1.2) and insist much more on the small divisors case, whose existence part is 
useful to us. To determine the role of n, we make the change of variable z = M T y in ( 1.1 ). 
One obtains that v(z) := um(Mz) solves 

-V • B(Mz)Vv = 0, z d > a 

v = v (Mz), z d = a 



(2.1) 



The family of matrices B = B a ^(y) £ Mn (M), indexed by 1 < a,/3 < d, satisfies for every 
i, j £{!,... N}, 

Bij = MA ij M T . 
and is hence 1-periodic, elliptic and smooth. 

2.1. Rational case. This case has been studied by many authors (see |AA99| IMV97j ). 
The assumption n £ QW 1 simplifies much the existence and convergence proof. Indeed, as 
n has rational coordinates, one can choose an orthogonal matrix M with columns in IRQ d 
sending e d on n. Subsequently, there exists a d-uplet of periods (L\, . . . L d ) such that 

z i-> B(Mz) and z v Q {Mz) 

are (L\, . . . L^)-periodic functions. Without loss of generality, let us assume that L\ = 
. . . = L d = 1. We can also fix a = for the moment. We come back later to this hypothesis. 
Then, lifting the Dirichlet data vo using ip £ C£° (M), compactly supported in [—1,1] such 
that < if < 1 and tp = 1 on [— h] , yields that v := v — ip(zd)vo(M.(z' , 0)) solves 

/ -V-B(Mz)Vv = V ■ B(Mz)V(^(z d )v (M(z',0))), z d > 
1 v =0, z d = Q ■ 



An appropriate framework to write a variational formulation for (2.2) is the completion of 
C oo (jd-1 x witn reS p ec t to the L 2 (T^ 1 x R + ) of the gradient. The fact that the 
source term V • B(M.z)'V (ip(zd)vo(M(z' , 0))) in (|2.2[) is compactly supported in the normal 



direction, allows to resort to a Poincare inequality. We can prove the existence of a weak 



solution v to (2.2) by the mean of the Lax-Milgram lemma. 

The existence part asserts that Vv £ I? (T d_1 x We aim at showing that 

actually decays faster. The key observation is that for any k £ N, := v(z',Zd — k) 
defined for Zd > k solves 

-V • B{Mz)Vv^ =0, z d > k 

=vo(Mz), z d = k ' 

Hence, one has for all k £ N the St-Venant estimate 

(2.3) \\^ v \\L 2 (T d - 1 x]k+l,oo[) - C [W^ V \\L 2 (T d - 1 x]k,oo[) ~ II (T^- 1 x]fc+l,oo[) J ' 

which yields the exponential decay of ||Vu || i2 ^T [J - 1 x]t oo[)' wnen * ~~ ^ 00 • The existence 
of the boundary layer tail £ M. N comes from the fact that fjd-i^ v (t)dt is a 

Cauchy sequence. The decay of higher order derivatives follows from elliptic regularity 
(see [ADN64J). This implies, through Sobolev injections, pointwise convergence of v(z', z d ) 
towards v,! , at an exponential rate. 

Let us come back to the assumption a = 0. Let a be any real number, v the associated 



solution of (2.1 ) and denote by a = a — [a] its fractional part. Then, v a = v(- + aed) 

-V • B(M(z + ae d )Vv a =0, z d > 

v a = v (M(z + ae d )) , z d = ' 



If we now assume that N = 1, d = 2 and that B(M-) is the constant identity matrix I2, 
we can carry out Fourier analysis to compute the tail. We get that the tail of v and 
v a is equal to 

-1 







vq (M(z' , a)) dz , 



which depends on a. Thus, the boundary layer tail is not independent from a in this rational 
setting. We now summarize the preceding results (forgetting about the assumptions L\ = 
. . . = Ld = 1 and a = 0) in the: 

Theorem 2.1 (lemma 4.4 in |AA99j . appendix 6 in |MV97] ). Assume that n G RQ d . 
Then, 

(1) there exists a solution v G C°° (IR ' -1 x [a, 00 Q of (2.1), unique under the condition 
that for all R> 

G L 2 ((-i2,fl) d - 1 x]a,cx3[), 

(2) anc? k > 0, v?j°° G depending on a such that for all a G N d ; /or a// 2 = (2/, 2^) G 
rod— 1 ■ 



"X Jo, 00 L 



(2.4) 



^-^l^^z)-^ 00 )! <c a . 



We conclude this section by a remark. It is concerned with the exponent k = n n in (2.4). 



In |NR00| . the case of a two-dimensional layered media £l n fi is considered. It is shown that 
exponential convergence of Vvu to 0, uniform in n G MQ rf , cannot be expected. Indeed, 
depending on n, K n can be arbitrarily small. 



2.2. Small divisors case. All the results we recall here stem from the original article 
|GVMll] by Gerard- Varet and Masmoudi. The periodic framework in the rational case 
makes the study of (2.1) more simple for the main reason that it yields compactness in 



the horizontal direction. One can thus rely on Poincare-Wirtinger inequalities, which are 



essential to prove the St-Venant estimate (2.3). In the non rational case, i.e. when n ^ 



one also attempts to recover a periodic setting. Note that for all z 



Iz'a) G 



pd-l 



x {a}, 



v (M(z', a)) = v (Nz' + M(0, a)), where N G M dA _ x (R) is the matrix of the d - 1 first 
columns of M. Therefore, Uo(M-) is quasiperiodic in z', i.e. there exists Vq = Vo(0,i) G M. N 



x {a}, v (M(z',a)) 



id-l 



defined for 6 G T and t g 1 such that for all z 

Vb(NV,a). Similarly, there exists B = B(6,t) such that for all z = (z',Zd) G 
B (M(z' , Zd)) = B(Nz' , Zd)- Hence, one looks for a quasiperiodic solution of (2.1); for 



details concerning quasiperiodic functions the reader is referred to [JK094J, section 7.1. 
Assume that there exists V = V(9, t) G M> N defined for 9 G T d and t > a such that for all 

z = (z',z d ) G R^x^oof, 



v(M(z',z d )) =V(Nz',z d ). 



Now, if V solves 



(2.5) 



N T V ( 
d t 



■B(6,t) 



N T V 
d t 



V 
V 



0, t > a 
Vq, t = a 



then v = v(z',Zd) = V(Nz',Zd) solves (2.1). 



We focus now on system (2.5). We examine successively the existence of a solution 



and the convergence towards a constant vector field. By making the change of unknown 
function, we have gained compactness in the horizontal direction, but lost ellipticity of 



the differential operator. System (2.5) is however well-posed (see |GVMll] proposition 2): 
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there exists a unique variational solution V G C°° (T d x]a, ooQ suc/i i/iai /or all / > 0, for 
all a G N d , 



(2.6) 



ft 



L 2 (T d x]a,oo[) 



J J°°(\N T Vod l t d$V\ 2 + \d\ +1 d^vf) dOdt < oo. 



To compensate for the lack of ellipticity, one uses a regularization method. We add a 
corrective term —tA to the operator. The regularized system is 



(2.7) 



^).B(fl,t)(^)^-^ =0, t>a 

V" =Vb, t = a 



The regularizing term yields ellipticity of the differential operator. Lifting the boundary 
data Vq into a function compactly supported in the direction normal to the boundary, 
makes it possible, as in the rational case, to prove the existence of a variational solution 



V L to (2.7). Carrying out energy estimates on (2.7) gives (2.6)-like a priori bounds on 
V L , uniform in i. Thanks to this compactness on the sequence (y t ) / >0 , o ne can extract 
a subsequence, which converges weakly to V variational solution of (2.5) when i — > 0. 



Uniqueness follows from the a priori bounds. An important point is that the justification 
of this well-posedness result does not resort to the small divisors assumption. 

We come to the asymptotical analysis far away from the boundary. Sobolev injections 
and the a priori bounds (2.6) yield pointwise convergence to of N T VgV(6, t), dtV(9, t) and 
their derivatives, when t — > oo, uniformly in 6 G T d . Let us investigate this convergence in 
more detail. Without loss of generality, we assume temporarily that a = 0. The idea is to 

^" V. In the same spirit as in the rational case, 



establish a St-Venant estimate on 
we look at the quantity 



d t 



K{T) := f f 
defined for T > 0. One proves that 

(2.8) K(T) < C (-K'(T)) 



N T V e V\ 2 + \d t V\ 2 )d9dt 



\V(0,T)\ d6 



with V(0,T) := V(8,T) — Jj d V(-,T). The stake is to control the second factor in the 
r.h.s. of (2.8). The key argument of [GVMllj is a type of Poincare-Wirtinger inequality 
implied by the small divisors assumption. Assume that n ^ RQ d satisfies (1.12). Then, 
for all 1 < p < oo, there exists C p > such that 



(2.9) 



\V(6,T)\ d6 < c p 



\N T v e v(e,T)\ do 



Note that the r.h.s. does not involve the L 2 (T d ) norm of V#U, but the norm of the 
incomplete gradient N^VgU. We give a sketch of the reasoning leading to (2.9). The 



proof relies on Fourier analysis in the tangential variable. In order to emphasize the role 
of the small divisors assumption, let us give an equivalent statement of (1.12): taking for 
the d - 1 first columns N G M^-i (R) of M, fll.l2| becomes 



»i 



(2.10) 



N J e \>c\i\ 



-d- 



Let T > be fixed, 1 < p < oo and p' its conjugate Holder exponent: i + h = 1. By 

Parceval's equality we can compute the I? (T rf ) norm on the l.h.s. of (2.9) and then apply 
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Holder's inequality 



with a 



\V(9,T)\ d6 < 



V?6Zd\{0} 




i B := \ and 7 := ^ 



,^ +T -, \ . For the first factor in the r.h.s., which 



p'(p-i)- 

represents the norm of V(-,T) in an homogeneous Sobolev space with negative exponent, 
we have, using (2.10): 



N T V e V(6,T) 



L 2 (T d ) 



2 bivry^r) |N T e| 2 >c \ V ^ T ) 



Sez"\{o} 



1 



2(d+T) ' 



This bound and (2.6) imply in particular a control on ||V(-, T) ||# s ^ T d^ f° r an s > 0. Hence, 

one bounds the second factor by a constant independent from T. 

The bound (2.10), in combination with (2.8), yields for every 1 < p < 00 the differential 

p+t — 

inequation on K, K(T) < C P {—K'(T)) 2 p from which we get the decay: for all T > 0, 

£+1 

< K(T) < C p T 1 ~p. We do the same on higher order derivatives considering for s G N, 



Ks(T) := 

0<\a\+l<s 



N T <9^V F| 2 + \d l t +1 d$V\ 2 ) dOdt. 



The existence of a boundary layer tail follows from a procedure similar to the rational case. 
For arbitrary a G M, one can prove (see proposition 5 in [GVMllJ) that the boundary 
layer tail does not depend on a. This fact is a characteristic of the non rational setting: 
Gerard- Varet and Masmoudi prove it under the small divisors assumption. Note that we 
manage to free ourselves from this assumption in section 6.2 To summarize (see |GVMll] 
propositions 4 and 5): if n ^ IRQ rf satisfies the small divisors assumption (2.10), then there 
exists a constant vector field G M. N , independent from a, such that for all m G N, for 
all a G N d , I G N, for all t > a, 



sup 

0eT d 



Q. + \t-a\ m )d?dl(V{6,t)-vR) 



< C„ 



We end this section by a translation of the existence and convergence statement for V 
into a statement for v solution of (2.1). We recall that for all z = (z',Zd) G M d_1 x]a,oo[, 
v{z',z<}) = V(Nz',Zd). A solution v of (2.1) (or vu of (1.1)) built like this is called a 
variational solution. 



Theorem 2.2 ( |GVMllj ). Assume that n £ 
(1) Then, there exists a solution v G C°° I 



(2.11a) 
(2.11b) 



00 



X 
t— >oo 



\ d z d v(-,t) 



[a, 00 [) of (2.1) satisfying 


dt < 00. 



(2) Moreover, if n satisfies the small divisors assumption (1.12), then for all m G N, 
for all a G N d , for all z = (z', z d ) eP"' 



a 



- x\a, oo[, 

\0? (v(z)-V%)\<C a , r , 



Estimates (2.11a) and (2.11b) rely on (2.6). Indeed, for every k G N, 

t— >oo 



l NTV ^||^(Tdx]t,oo[) + \\ d t V \\m(Td x] t,oo[) 
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which together with Sobolev's injection theorem yields (2.11a) for k sufficiently big. As 

) Lf(]0,oo[) 



\d t V\\ Hk 



H k (T d x]a,oo[) 



> 



a similar reasoning leads to (2.11b). The a priori bounds (2.6), contain further information: 



for k > sufficiently large and k' > 1, 



(2.12) 



0, 



V 



> 



H k '+ k - 1 (T d x]a,oo[) 



d. 



V 



H h g^ d ) 



> c 



dt 



V 



L e°( Td ) 



> c 



V k 'v 



L?(]a,oo[) 



Ll d {]aM) 



We resort to the latter in section l5~2l with k! = 1 or 2. 

3. Poisson's integral representation of the variational solution 
One can always assume that a = 0. We do so for the rest of the paper (except in section 



6.2), even if it means to work with := Vbi(- + an) instead of vu solving (1.1). The 
main advantage of this assumption is that the domain £l n = {y ■ n > 0} is invariant under 
the scaling y i— )■ ey for e > 0. The purpose of this section (see in particular section 
is to prove uniqueness for the solution to (1.1) in a larger class than that of jGVIV 



3.2) 



HI 



This result (theorem 3.5 below) allows to represent the variational solution by the mean 
of Poisson's kernel. 

3.1. Estimates on Green's and Poisson's kernels. Let G = G (y,y) £ Mjv(M) solving, 
for all y G f2 n , the elliptic system 



(3.1) 



-V y -A(y)V y G(y,y) 



%-y)Iiv, y-n>o 

0, y ■ n = 



I G(y,y) 

with source term 5(- — y) (5(-) is the Dirac distribution). The matrix G is called the Green 
kernel associated to the operator —V • A(y)X7- and to the domain £l n . Its existence in f2 n 
is proved for d > 3 in [HK07J (see theorem 4.1), and for d = 2 in [DK09| (see theorem 
2.12). Similarly, G* = G* (y,y) G M/v(M) is the Green kernel associated to the transposed 
operator —V • A*(y)X7- and the domain f2 n , A* being the transpose of the tensor A defined 
for all a, f3 £ {1, . . . d} and for all i, j G {1, . . . N} by 



(A*) 



a/3 



A 



/3a 



(A /3a ) T 

From the smoothness of A and local regularity estimates, it follows that 

G G C°° (fl n xn n \{y = y}) . 

Moreover, the following symmetry property holds: for all y, y G f2 n , 

G T (y,y) = G*(y,y). 

Using Green's kernel, one defines another function, the Poisson kernel P 
M N (R): for all i, j G {1, . . . N}, for all y G Q n , y G 3ft n , 

(3.2a) Pij(y,y) ■= - A tj(v) d y a G ik{y,y)np 

(3.2b) =-(A*)%(y)a Sa G* ki (y,y)n f) 

(3.2c) =-\A^ a (y)dy a G*(y,y)n p 



P(x,x) G 



(3.2d) 



[^(y)V S G*(y,y)-n]i. 

12 



The kernel P* = P* (y, y) G Mn(R) is denned in the same way as G* . 

If one considers Green's kernel G° = G°(y,y) G Mjy (M) associated to the constant 
coefficients elliptic operator —V • ^4°V- and to the domain f2 n , there exists C > (see 
[ADN64J section 6 for a statement and |SW77| V.4.2 (Satz 3) for a proof) such that for all 
Ai, A 2 G N d , for all y, y G Q n , y ^ y, 

(3.3a) \G°(y,y)\ < C(|ln|y - y|| + 1) , if d = 2, 

(3.3b) |G°(y,y)|< , C |d _ 2 , if d>3, 

y - y 



(3.3c) d^d^G°(y,y) 



^ r^W ' if|Aii + |A 2 |>i. 



One has similar estimates on Poisson's kernel P° = P°(y, y) G M/v (ffi) and its derivatives. 
Such bounds on Green's kernel and its derivatives are not known for operators with non 
constant coefficients. Let us temporarily assume that neither (A2) nor (A3) hold. We 
know then from |KK 10| (see theorem 3.3) that under certain smoothness assumptions 
on the coefficients there exists R ma x G]0, oo] and C > 0, such that for all y, y G O n , 
|y - y\ <Rmax implies 

(3-4) \G(y,y)\< . 

\y - y\ 



If Rmax < oo, estimate (3.4) does not provide a control of G(y, y) for y and y far from each 
other. However, under the periodicity assumption (A2), we have the following improved 
global bounds: 

Lemma 3.1 QAL87] theorem 13 and lemma 21, |GVM| lemma 6). There exists C > 0, 
such that 

(1) for all y, y G Q n , V / V, 

3.5a) \G(y,y)\ < C(|ln|y - y|| + 1) , if d = 2, 

3.5b) \G(y,y)\<- ifd>3, 

\y - y\ 

and for all d > 2, 

3.50 \G(y,m<C b '- n) % n \ 

\y - y\ 

3-5d) \V y G(y,y)\ < C -^, 

\y-y\ 

3.5e) | V , G(y , i> )|< C f^ t+ to-"«^A 

\\y-y\ \y-y\ J 

3.5f) | V ,V^G(y,y)|<— 

\y - y\ 

(2) for all y G £l n , for all y G dVt n , for all d > 2, 
3.6a) | P(y ,y)|<C7^^, 

\y - y\ 

3 - 6b ) \Vy p {yi y)\ < c \-^— d + , y 'Z +1 ) ■ 

\\y-y\ \y-y\ J 

By continuity of G and its derivatives, up to the boundary, the estimates on Green's 
kernel naturally extend to y ^ y G £l n . These bounds are of constant use in our work: for 
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a proof in the half-space see |GVM| (appendix A). The key argument is due to Avellaneda 
and Lin. The large scale description of G boils down to an homogenization problem. In 
the paper |AL87j . the authors face this homogenization problem under the periodicity 
assumption (A2) and manage to get uniform local estimates on u £ = u £ {x) satisfying 



(3.7) 



-V- A(f) Vu £ 

u £ 



f, 13(0,1) 

g , xer(o,i) 



We recall the two local estimates useful in the sequel: 

Theorem 3.2 (local boundary estimate, |AL87j lemma 12). Let fj,, 5 be positive real 
numbers, fj, < 1, F G L d+S (D(0, 1)) and g G C ' 1 (T(0, 1)). Assume that f = V • F. 
There exists C > such that for all e > 0, if u £ belongs to L? (D(0, 1)) and satisfies (3.7), 
then 



(3.8) 



KHco>m(d(0,!)) ^ c 



lL 2 



(o(o,i)) + \\ F \\L d + s (D(o,i)) + ll5 f llc°> 1 (r(o,i)) 



Theorem 3.3 (local boundary gradient estimate, |AL87j lemma 20). Let u, 5 be positive 
real numbers, v < 1, f € L d+S (D(0, 1)) and g G C 1 '" (r(0, 1)). 



There exists C > suc/i that for all e > 0, i/n e belongs to L°° (D(0, 1)) and satisfies (3.7), 



(3-9) l|Vul Loo(D(0i|)) < C 



u IIl°°(D(0,1)) + ll/lli d + 5 (£»(0,l)) + llsllc 1 '"^^,!)) 



3.2. Integral representation formula. We aim at showing that the solution vu of (1.1 ) 



can be represented in terms of an integral formula involving Poisson's kernel. One of the 



main difficulties arises from the fact that the Dirichlet data vq of (1.1) is not compactly 



supported in the boundary. The function v = v(z) G K , such that for all z G 



pd-l 



v(z) '.= Vf,i(M.z), solves (2.1). Let us now recall what we consider as a solution of (2.1). In 



fact, we have two kinds of solutions. The first corresponds to the variational construction in 



[GVMllJ: see section 2.2 in particular theorem 2.2 Poisson's kernel P = P{y, y) G Mjv(R) 
associated to the domain Vt n and the operator —V • A(?/)V- makes it possible to define a 



second function wu = Wbl{y) £ K solving (1.1). For all y G Cl n , we define Wbi(y) by 



wu(y) = P(y,y)v {y)dy 

' y-n=0 



and w by for all z G M. d 1 x]0, oo[ 



w(z) := w b i{Mz) 



P(Mz,y)v (y)dy = / P(Mz, Mz)v (Mz)dz. 



Proposition 3.4. The function w belongs to C°° (R 1 x [0, oo[) and satisfies (2.1). Fur- 
thermore, 



(3.10a) 
(3.10b) 



\V z w\ 



t—>oo 



lL°°(R d - 1 x]t,oo[) ' 
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Let us give a sketch of how to deduce these properties from the bound on V y P given in 



-x 0,oo , 



fl3.6[ ). For all z € K d_1 ; 

\V z w(z)\ < / |M T Vj / P(Mz,Mi)f (Mi)| dz 

Jz d =0 



< C 



< C 



< C- 



Lz — z\ d 



+ 



1 



dz 



+ 



,z 



d + \z'-~z'\ 2 Y + 

1 1 

7 + 



d+1 



d+1 



dz 



dz 



*d jRd-l ^(1 + | Z / _ z'\ 2 )2 (1 + | z ' - Z'\ 2 )~ 



from which one gets (3.10a) as well as (3.10b). 

Our goal is now to show that the variational solution and the Poisson solution coincide. 

Theorem 3.5. We have v = w. 

We work on the difference u := v — w, which is a C°° solution of 



oo {md-l. 



- V • B(Mz)X7u =0, z d > 
u =0, z d = 

We intend to show that u is zero proceeding by duality. Let / G C\ 
take U = U(z) S the solution to the elliptic boundary value problem 

-V-B*(Mz)VU = /, z d >0 
U =0, z d = 

given by Green's representation formula: for all z £ W 1 ^ 1 x]0, oo[, 



]0,oo[). We 



I7(z) 



G*(Mz,Mz)f(z)dz. 



z d >0 



The bounds on G* (3.5c) and its first-order derivative ( |3.5d ), yield: 

Lemma 3.6. There is C > such that for sufficiently large z € M.^ 1 x]0, oo[, i.e. far 
enough from the support of f, 

Zd 



(3.11a) 
(3.11b) 



\U(z)\<C- 



\VU(z)\ < C- 



zl + \z'\ 2 ) r 



\zi + \z 



x d ' 

l\2\ 2 



Moreover, thanks to (2.11b) and ( 3.10b[ ), one manages to estimate u in L°°: there is 
C > 0, such that for all zir -1 x]0, oo[, 



2d 



d z .u{z',t)\dt<z} 



(3.12) \u(z)\ < 

We now carry out integrations by parts: 

u(z)f{z)dz = - 



d Zd u(z,t) dt) <Cz> 



z d >0 



u{z)V ■ B*(Mz)VU(z)dz 



z d >0 



B a/3 (Mz)d Z0 u(z)d Za U(z)dz 



z d >0 



z d >0 



8 Za (B<*P{Mz)d Zfj u(zj) U(z)dz = 0. 
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To fully justify the preceding equalities we have to integrate by parts on the bounded 
domain [-R^]^ 1 x [0,R] and prove that the boundary integrals vanish in the limit 
R — > oo. We actually show that 



(3.13a) 
(3.13b) 



d([-R,R] d - l x[Q,R]) 
f 

9([-ii,il]< i - 1 x[0 ) ii]) 



u (z) (B* (Mz) VU (z)) ■ n(z)dz 0, 



[B (Mz) Vu (z)] ■ n(z)U (z) dz 0. 



On the one hand, (3.12) together with (3.11b) yields 



a([-il,fl] d - 1 x[0,.R]) 

< c 



u (z) (B* (Mz) VU (z)) ■ n(z)dz 



\ a~ ^ n T>d-l 



d I „|<i 



d([-R,R] d - 1 x[0,R])\([-R,R] d - 1 x{0}) \ z 



dz < C r R 

R d ~2 



c_ 

r > 

R2 



which gives (3.13a). On the other hand, it follows from (2.11a), (3.10a) and (3.11a) that 



R^^xiR} 



[B (Mz)Vu(z)]-n(z)U(z)dz 



< C 

< C 



R,R] d ~ 



\Vu(z',R)\ \U(z',R)\dz' 



iVnl 



R 



- 1 x]R,oo[) 



(R 2 + \z'\ 2 ) 



-dz' 



< C ||Vu|| ioo j- Rd _ lx]i j i0o Q — > 



and that 



{R}x[-R,R} d - 2 x[0,R] 



[B (Mz) Vu(z)]-n(z)U(z)dz 



< C 

< C 



IVnl 



-R,R] d ~ 2 x[0,R] 



Wu\ 



-l,l] d - 2 x[0,l] 



L x]2; d ,oo[ 



- 1 x]Rz d ,oo[) 



(R? + N 2 + • • • M 2 ) 

(l + \z 2 \ 2 + ... \z d \ 



-dzo 



d """2 

2 



• dz d 
dz 2 . ■ ■ dz d 



tend to by dominated convergence, which yields (3.13b) and terminates the proof of 
theorem 13.51 



For practical convenience, we have argued that v = w. Yet, theorem 3.5 proves that 
the variational solution v^i of (1.1) equals the Poisson solution Wyi- This allows to work, 
for the rest of the paper, with the solution v^i of (1.1) satisfying (1.14a) and (1.14b), no 
matter whether this solution is constructed variationally or via Poisson's kernel. Thanks 
to the bound (3.6a), Vf>i is seen to be bounded on O n . 

3.3. An homogenization problem. Studying the tail of i>bi, i.e. the limit when y-n —> oo 
of Vbi(y), boils down to describing the asymptotics of P(y,y) for y far away from the 
boundary dVt n . One of the main focus of our paper is thus to expand P(y, y) for \y— y\ 3> 1, 
where y S Vt n and y £ dU n , i.e. to describe the large scale asymptotics of P. Let y E fi n 
and y £ dQ n and 

1 

£ := ; r. 



\y - y\ 
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If \y — y\ 3> 1, then e <C 1 is a small parameter. Introducing the rescaled variables 

x := ey G Cl n and x := ey G Q, n 

yields \x— x\ = 1. Such a scaling transforms our initial question of the large scale asymptotic 
description of P into the study of -^rP (f > § ) f° r £ ~ * and |x — x\ close to 1. 

Lemma 3.7. Let e > and ca// G e (resp. P £ ) the Green (resp. Poisson) kernel associated 

to the operator L e = — V • vl (|) V- and the domain Q n . 

Then, 

(1) /or a// x, x G fi nj 



(3.14) 



G £ (x,x) 



r d-2 



G 



' / 1 

e e 



(2) for all x G £l n , x G 



(3.15) 



P e (x,x) 



1 / r ,r 



e d - 1 \£'e 



Proof. This lemma follows easily from Green's integral representation formula. Let / 6 
C~ (n n ), u e = u e (x) G 1^ the solution of 

' -V-A(f)Vn £ = £/(§), x-n>0 
u e = 0, x ■ n = 



and it = ii(y) G R the solution of 



-V • A(y)Vu = /, y • n > 
u =0, y • n = 



For all x G Q n , 



G £ (x, x)-^f \^-jdx = u £ (x) = u ^— ^ 

G^y)f{y)dy 



— G[ - - 



dx. 



which yields (3.14); (3.15) easily follows from analogous ideas. 



□ 



It immediately follows from the definition of G e , that for all x G Q n , G £ (•, x) solves the 
system 

-V x ■ A (|) V X G £ (x,x) =5(x-x)l N , x-n>0 
G £ (x, x) = 0, x ■ n = 

The Poisson kernel P e satisfies: for all i, j G {1, . . . N}, for all x G f2 n , x G 9f2 n , 



(3.16) 



(3.17a) 
(3.17b) 
(3.17c) 

(3.17d) 



P! j (x,x) = -A«h^\d Sa GUx,x)n f3 



(Anin^d^Gif^x)^ 



x 



1 ( ~ Wg^x)^ 
^ ( | ) V x G*> £ (x,x)-n 



j i 



where G*' £ is the Green kernel associated to the operator L* ,£ = — V - A* (|) V- and the 
domain £l n . 
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The estimates (3.5) and (3.6) can be rescaled. In particular, there exists C > such 



that for all e > 0, for all x, x G £l n , i /i, for all d > 2 
(3.18a) 

|G £ (x,£)| < 



(3.18b) 
(3.18c) 
(3.18d) 



G £ (x,x)\ < C7(|ln|x-x|| + 1) , if d = 2 (see [AL87J theorem 13 (ii)), 
( if d > 3, 



I I 2 



|G £ (x,x)| < C 
|V £ G E (x,x)l < 



'x ■ n) (x ■ n) 

i ~ \ri : 

\X — X\ a 

C 



~ \rl— 1 ' 

x — x\ a 1 



and for all x G Q n , x £ d£l n , 



(3.18e) 
(3.18f) 



\P e (x,x)\ < C 
\V x P £ (x,x)\ < C 



x ■ n 



x — x\ a 



\x — x\ d 



+ 



x ■ n 



\x — X 



d+1 



According to lemma 3.7, we now deal with highly oscillating kernels G £ and P £ , instead 
of looking at G and P. Hence the asymptotic description of G (resp. P) at large scales, 
is replaced by an homogenization problem on G £ (resp. P £ ). This fact, which has been 
stressed by Avellaneda and Lin (see |AL91| corollary on p. 903), is the cornerstone of our 
method. 



4. Homogenization in the half-space 
This section is concerned with the asymptotics, for small e, of u £ = u £ (x) £ R N solving 



(4.1) 



-V-A(f)Vn £ =/, x-n>0 
u £ =0, x ■ n = 



The study of this dual homogenization problem is preparatory to the expansion of the 
Green and Poisson kernels in the next section. In the introduction, we defined the interior 
and boundary layer correctors to u £ up to the order e and reviewed some error estimates in 
the case of a bounded domain Vt. The purpose of this section is to show similar estimates, 
yet in the unbounded domain VL n . 

Let / G (fin)] note that the support of / may intersect the boundary. We define 
u° = u°(x) G ~R N as the solution of 



(4.2) 



-V • ,4 Vu = /, x ■ n > 
u° =0, x • n = 



u 1 = v}(x,y) G 1^ by u 1 (x,y) := X a (y)dx a u°(x), for all x G f2 n and y G T d , and 
u,f = u b f(x) G M. N as the Poisson solution to 

-V ■ A(f)Vulf =0, x-n>0 



uif — - a 



1 (x, f) = -X a (f) d Xa u°(x), x-n = 



The variational solution u £ (resp. u°) coincides with the solution given by Green's integral 
formula. Besides, u £ , u , as well as u b f belong to C°° (Cl n ), thanks to the smoothness of 
the boundary d£l n , using local regularity estimates from |ADN 64j. The rest of this section 
is devoted to the proof of the following proposition: 
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Proposition 4.1. Let r^f := u £ (x) — u°(x) — ex a (§) d Xa u° — eulf(x), and 5 > 0. 
Then, there exists a constant C > 0, such that for all e > 0, 

(4.3a) 

(4.3b) 
(4.3c) 

The proof of (4.3a) relies on estimates in L°° of r^' e solution of the following elliptic 
system 

-V-A(f)Vr^f = f £ , x-n>0 
0, x ■ n = 





l,s 
r 6i 




<Ce\\ 






l,s 


L°°(n n ) 


<cy 




u £ 


-«°| 




<Ce\\ 





' r lf 

Le 



where / e = f £ (x) £ is given by 



/ £ := / + V • A 



Vu° + eV • i4 



v X 



d Xa u°{x) 



The idea is to use the integral representation provided by Green's formula in order to 
bound r^' £ . However, as such, f e does not seem to be of order e. Let us thus work on f e 
to make its structure more explicit. Expanding f £ , we get for all x G r^ n 



(4.4) f(x) 



[V y ■ A(y)V x v° + V y ■ A^VyU 1 } 
[/ + V s ■ A(y)V x u° + V x • A(y)V y u 1 + V y • A{y)V x u l ] 
+ e[V x - A{y)V x v}-] 



We aim to get rid of terms of order e 1 and e° in (4.4). The term of order s 1 easily 
cancels thanks to (1.4): 

V, • A(y)V x u° + V, • A{y)V y u l = [d Va (A^(y)) + d Va (a^ {y)d y[j x i (y))] d Xl u° = 0. 

For the term of order e° in the r.h.s. of (4.4) 



(4.5) [/ + V x ■ A(y)V x u° + V* • A{y)V y u l + V y ■ A{y)V x u l ] (x, 

= [f + V x • v + V„ • A{y)V x u l ] 
where v = v(x, y) := A(y)V x u° + A(y)V y u 1 = A(y) [V x u° + V^u 1 ] , we notice that 



V y ■ (v - A°Vu ^ 



and 

As v — A°Vu° factors into $Vn°, with 



/ (v - A°Vu°) = 0. 



D N 2 xd 2 



$ = : = A(y) (I +V vX (y)) - [ A(y) (I +V yX (y)) dy G 
one can take advantage of the classical lemma: 

Lemma 4.2. There exists a smooth function ^> = ^(y) € R^ 2 ^ 3 such that for all y £ T rf 7 

$(y) = V y • *(!/). 

Via lemma 4.2 (4.5) becomes 

/ + V x • w + V y • ^(y)V x u 1 = V x ■ [v - A°V x u } + V y • A(y)V x u l 

= V x ■ (V y • (*(y)) Vu°) + V a • A(y)V :E n 1 
= V„ • (¥(y)V 2 u°) + ^(y)V 3 u°. 
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Subsequently, for all x G fi n , 

= [V, • (*(y)V 2 u )] 
= [V, • W)] (*, 



+ eV ■ [A 



VV) 



+ E [v I -(f(i,)vV)] 



.r 



+ eV • (A 



eV- * 



X 



.r 



x, - 

e 



x 



V V ) +eV- (A 



X 



V 2 u° 



VV. 



Hence f £ = eV • /i £ + eg £ , with 



.r 



/i £ = h £ {x) := * ( - ) V 2 u° +A 



vV, 



<f = <f(x) := 



vV. 



The next lemma contains estimates on h £ and g £ for large x. 

Lemma 4.3. There is a constant C > 0, such that for all x sufficiently large, for all e > 0, 



(4.7a) 
(4.7b) 



\h £ (x)\ < C 



1 



< C- 



1 



\d+l ' 



Proof. Let A G N d , 2 < |A| < 3, and R > such that the support of / is included in 
B(0,R). It follows from (|3.3[) that for x large enough, 



9 A n°(^)| < / \d£G"(x,x)\ \f(x)\dx 



Ar<0/ 



< c 

< c 

< 



1 



B{0,R) \X — X 



id-2+|A| 
1 



dx 



B(0,R) (\x\ - R) 

c 



d-2+|A| 



(\x\-R) 



d— 2+|A| 



o 



dx 



1 



\ x \d-2+\A\ 



which ends the proof. 



□ 



These preliminaries being done, we now turn to the estimation of r^' £ . Let x G Q n be 
fixed. Green's formula yields 



(4.8) 



rlf( x ) = e G £ (x, x) (V • h £ + g £ ) (x)dx. 

J Q n 



We concentrate on each term of the r.h.s. of (4.8) separately. The strategy in both cases 
is to split the integral in two parts: 

(1) for x close to x, one relies on Young inequalities to bound this part in L°°; 

(2) for x, far from x, one uses (4.7a) and (4.7b) to show that this part can be made 
arbitrarily small uniformly in x. 

Let R > and assume for the moment d > 3. An integration by parts 



G £ (x,x)V ■ h £ (x)dx 



{VxG £ ){x,x)h £ {x)dx, 
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together with (3.18b) and (3.18d) gives 
r\f (x)\ < e f ° \h £ (x)\ dx + e I " \g £ (x)\ dx 



< £ 



n„ \x - x 
C 



C_ 

n n \x-x\ 



+ £ 



\x — x\ d 
C 



i^-B(o,R)( x ~ x)\h e (x)\dx + 



C 



\x — x\ 



^B(o,R)4 x ~ x ) \h £ (x)\dx 



x — X 



c 



(x - x) \g £ {x)\ dx + e I - 1 B (o,r)4 x ~ x ) \F(x)\ dx, 

J R d \x - x\ a 

where h e (resp. g £ ) is the extension of h £ (resp. g £ ) to M. d by outside of Q n . Let us first 
concentrate on the terms involving h £ . First, it simply follows from lemma 4.3 that |/i £ (x) 

is O ( i^jdj i n a neighbourhood of oo. One can find p, p! > 1 such that 



1 1 d 

p p a — 1 



and p > 1. 



Therefore, 



1 



x 



T 1 B(0,H) c (^) 



R—too 







LP( 



and h £ is bounded uniformly in e in L p ' (R d ) . Thanks to Young's inequality, 
C 



(4.9) 



/ 



x — X 



Zi^-B(o,R) c (x - x)\h £ (x)\dx 



< 



\h £ 



\Lp' 



uniformly in x 6 O n . The r.h.s. of (4.9) can be made less than 



R— >oo 



for an R > 



large enough. We now carry out the analysis of the integral on \x — x\ < R. An adequate 
choice of the exponents in Young's inequality leads to (4.3a). Indeed, for all 1 < q < -t^j. 



n~p=i ^B(o,R) ^ L q . From the condition 1 = | + on Young exponents, one deduces 
q' > d. Yet, for all 5 > 0, for all q' := d + 5, thanks to elliptic regularity and Sobolev's 
injection 



If 



\Li'(n n ) 



< C||V 2 u || r „,, o , < C||VV|| . d+s 

on 



\Li'(n n ) 



<C\\u 



< C 



Young's inequality finally gives 



1 



\x — x\ d ' 



zi l B(o,R){x - x)\h £ (x)\dx < C 



The reasoning for g £ is almost the same, except for the exponents in Young's inequalities 
which have to be adapted. 

Let us briefly indicate how to treat the case d = 2. Each term can be estimated as 
above, except the ones involving g £ . As before, we split the integral: 



G £ {x,x)g £ {x)dx = I G £ (x,x)l D(0:R )(x - x)g £ {x)dx 



+ G £ (x,x)l D ( 0tR )c(x - x)g £ (x)dx. 
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For x close to x we bound G £ by (3.18a 



G E (x,x)l D , 0R \(x - x)g £ (x)dx <C (|ln|x - x\\ + 1) 1 B (o,r)(x - x)g £ (x)dx, 
and for \x — x\ > R we have recourse to (3.18c) instead of ( |3.18b ) 



G £ (x, x)l D ^ R) c(x - x)g £ (x)dx 



< C 



\x • n\\x • n\ 



\x — x\ 2 



1 D(o,R)4 x ~ x ) \9 £ (x)\ dx 



< C 



\(x — x) ■ n\\x ■ n\ \x ■ n 



\x — x\ 2 



+ 



x — x\ 2 



1 B(o,R)4 x ~ x ) \9 £ (x)\ dx 



< C 



\x — X 



^B^R^ix - x)\x\ \g £ (x)\ dx 



+ 



;1b(o,r)4 x ~ x )\ x \ 2 \a £ (x)\ dx 



\x — x\ 



and estimate these terms, uniformly in x and £, via Young's inequality. The bound (4.3a) 



on r^' £ is shown. 



As elliptic regularity and Sobolev injections imply 



<C Vu 



o 



< C 



w 2 'i +s (n n ) ~ " UJ l V 1 >i+ a (n») ' 



it remains to establish (4.3b) in order to prove (4.3c). Let x G Q n . Poisson's representation 



formula for uh £ yields 



u l b f{x) 



[ P%x,x)x a (-)d Xa u u (x)dx. 



From the bound (3.18e) on P £ , one gets 
\u b f(x)\ < [ \P £ (x,x) 



dx 



< C 

< C 

< c 



dQ r , 



x ■ n 



\x — x\ 



1 fx 2 



+ \x' — x'| 2 ) 



1 



1 i l~/|2\2 

1 + \x n ' 



<C Vu 



< C Vu° 



o 



| d Xa u° (x) | dx 

£d ] d Xa u° (M (x', 0)) I dx' 

d Xa u° (M (x' -x d x',0))\dx' 

-dx 



1 (1 + |X'| 2 )5 



< C 



with x := M x. This establishes (4.3a) and proposition 



4.1 



5. Asymptotic expansion of Poisson's kernel 



We intend to get the asymptotics of P £ = P £ (x,x), defined by (3.17), for x in a neigh- 
bourhood of the boundary d£l n and \x — x\ close to 1. Our method is in three steps: 

(1) homogenization of u £ solution of (4.1) for / G (f^ n ); 

(2) expansion of G £ , thanks to a duality argument; 

(3) expansion of P £ via (3.17) and the expansion for G £ . 

The first point has been the purpose of section [4] We now turn to the second point. 
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5.1. Back to Green's kernel. 



Proposition 5.1. For all < k < ^, there exists C K > 0, such that for all e > 0, for all 

x, x G O n , \ < |x — a; | < 4 implies 

\G £ (x,x) - G°(x,x)\ < C K e K . 



We start from (4.3c) and proceed by a duality method, mimicked from |AL91| . to prove 
the estimate on the kernels. The key is as usual Green's representation formula. Let 
x G £l n and e > be fixed for the rest of the proof. We look at 

o £jX := sup \G £ (x,x) - G°(x,x)\ . 

f<|z-x|<4 

The l.u.b. o £)X is reached for at least one x E)X , which may be on the boundary dtt n . From 
(3.18d), one obtains the existence of C\ > 0, independent from e and x, such that for all 

X £ ^tyij g | X X^ ^ 5j 

\V s G £ (x,x)\ + \V x G°(x,x)\ < Ci. 
Let p S:X := 2N^Ci ' ^ ne can a l wavs increase Ci, so that p EjX < 1 and 

1' 



B (x £ , x ,Pe,x) C B(x, 5)\B [ x, 



Then: 



Lemma 5.2. There exists i, j G {1, . . . N} such that for all x G D (x e ,x, Pe,x), 

0~E,X 



\G! j (x,x)-G^(x,x)\ > 



2N 2 ' 



Proof. From x £>x ) — G°(x,x £iX )\ = o £:X it comes the existence of i, j G {1, . . . N} 

such that 



\Gij(x,X £jX ) G i j(x, X £) x)\ — jy 2 " 



(5.1) 

The integers i, j are now fixed such as ( |5.1| ) holds. For all x G D (x £iX , p £tX ), either 
GUx,x)-G%{x,x) 



< 



> 



27V 5 : 



or 



Gf,(x,x)-G?,(x,x) 



< glrf' Assume the latter. Then, 



< \Gfj(x, x £)X ) — G®j(x, x £)X )\ — \Gfj(x,x) — G®j(x,x)\ 



— \Gij(x,x £)X ) — G 6 j(x,x)\ + \G®j(x,x £ 



Gij(x, x)\ 



< 



\V X G £ 



+ N X G L 



lL°°(i<|a;-x|<5) T II v x ^ \\ L°° (\<\x-x\<S) 



< C\p £jX 



2N 27 



and finally 



\GUx,x)-GUx,x)\ > 



Or 



2N 2 ' 



□ 



Take ip G (B(x £)X , 1)), with < ip < 1 and ip = 1 on B (x £jX , |). Note that the 
support of p may intersect the boundary <9$7 n . For p > 0, we define ip p by ip p := cp G 
C~ (5 (x £ja; , p)); we have: 

C 



< C and 



PllL°°(fi„) 



< 



For p = p £jX , the constants above do not depend on e or x. 
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Let i, j the integers given by lemma 5.2 The intermediate value theorem implies that 
Gfj(x, x) — G®j(x, x) has a constant sign on D (x £ ^ x , p £jX ). Up to the change of / in — / in 
what follows, one can always assume that GSjfax) - G^(x,x) > 0, which automatically 
yields 

G?.(x,x)-G°.(x,x)> a 



'e,x 

2iV2 



for all x G D {x £)X , p £ ,x)- We now carry out the duality argument. For this purpose, take 
/ = f(x) := ip Pl , x (x)ej, where ej is the j th vector of the canonical basis of l w and denote 
by u £ = u e (x) G R N and u° = u°(x) G R N the solutions of ( |4~Tj ) and ( |4~2| ) with r.h.s. 
/. We remind that /, u £ as well as u° depend on x. Thanks to Green's representation 
formula, 



[u%x)-u°(x)l 



L 



[(G%x,x)-G°(x,x))f(x)] i dx 

[Gtj(x,x) - Gfj(x,x)] <p Pex (x)dx 



> 



> 



Gfj(x,x^) — G^j(x,x^dx 



a F 



Yet, we know from (4.3c|, an estimate of if — u°: 
.on 



dx > Cp 



d+i 

£,X ■ 



C s e 



d-2S 

<C s ep d £ % 2S . 



Putting together these bounds, we get 



Cpi% 1 < \[u%x)-u°(x)] i \ < \u%x)-u°(x)\ < \\u £ -u \\ Lx{nn) < C 5 e Pi 



d-2S 
d+26 
S.X > 



which summarizes in 

£ > ^SPe,x = (^SPe,x 

for every 5 > 0, with a constant C$ independent from e and x. The inequalities 



(T £ ,x < C 5 p £ , x < C s e d+ d+i 



contain the asymptotic expansion of G e at zeroth order of proposition 5.1 One can follow 
the same reasoning as above, changing A in A*, to obtain for all < k < -j, for all 

x 2 x £ ■) ~a_ \x X | 4^ 



(5.2) 



G*' £ (x,x) -G*'°(x,x)\ < C K e K . 



5.2. Homogenization of Poisson's kernel. Let < e < 1 and x G O n be fixed. Assume 
that x is close to the boundary, say x ■ n < 4 to fix the ideas. According to ( |3.16 ), G*' £ (-, x) 
satisfies 

-V x -A*(l)V x G*> £ (x,x) = 0, x€D(x,A)\D(x,\) 
G*' £ (x,x) = 0, x, G r(x, 4) \ r (x, j) 

This leads to the idea that one can apply theorem |3.3| to an expansion of G*' e for which 
a local estimate in L°° is known. Doing so, one has to handle carefully with the trace 
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on T(x,4) \ T (x, I). Take for example Z*> e > x = Z*> e > x (x) G M N (R) defined for all x G 
D{xA)\D(x,\) by 



(x) : = G*' e (x, x) - G*'°(x, x) - e X 



d Xa G*'°(x,x 



dx a dxpG*'°(x,x). 



Then, for all < v < 1, 

\\£*,S,X\ 



X 



dx a G*>°(x,x) + er 



,a/3 



d Xa d X/3 G*'°(x,x) 



(s-») 



ci."(r(*,4)\r(*,i)) 

which worsens our estimates. One way of getting around this difficulty is again to introduce 
a boundary layer term in the expansion . This term has to cancel the trace on the boundary 
due to the first-order term. 

For all 7 G {1, . . . d}, let G^' 1 ' 7 = G^' 1 ' 1 '(y) G M N (R) be the solution of 



-Vy-A*(y)VyG* b ^ = 0, yeU n 

-x*'Hv), ytdtt n 



in t he sense of theorem 2.2 or proposition 3.4, both being identical according to theorem 



3.5 



by 



We then define the boundary layer Gtl ' e = G^ ,1,£ : (x,x) G M/v(K) at first order in e 



G b i ,e (x, x) — Gf 



,1,7 / 



a £ 7 g*< u (x,x) 



for all x G O n , i/i Instead of Z*' £ ' x one focuses now on Z,f' = Z^ e ' x (x) G Mtv(K) 
such that 

/7*,£,2! 
■^1 



w (x) : =G*' £ (x,x)-G*' u (x,x)-e X * 



^ Q G*'°(x,x)- e G*^ £ (x,x) 



^^G*'°(x,x)-eV' a (-)G; 



-,*,!,£ / X 



d Xa d X0 G*>°(x,x 



for all x G Z)(x,4) \ Z) (x, |). The method is now similar to the one, which led to the 



estimate on r^' e : Ztf' x satisfies 



-V- A* (?) VZ,*' e 



• ~bl 



+ xGD(x,4)\ J D(x,|) 

_ £ 2 T *, a p ^ Q ^G*'°(x, x), x g r(x,4) \ r (x, I) 



where 



F £ := V £ • A* - V £ G*' u (x, x) + eV 2 • A* - V £ x 



x 



<9x a G*'°(x, x 



+ e 2 V x ■ A* 



Vx r 



d Xa d Xg G*'° (x, x 



and 



(5.3) ^:= £ V 2 -A*(^)V £ G^' £ (x) 



+ e 2 V £ • A* ( ? ) V, ( X *' a (-) G*^ ( ~- ) d^G*>°(x,x) ) . 
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Applying (3.9), one obtains 



(5.4) \\VZ-, 



bl \\l°°(D(x,3)\D(x,±)) - U 



z: 



l>l \\L™(D{x,A)\D{x,\)) 
+ \\ FS \\L d + s (D(xA)\D(x,\)) + W F bl\\L d + s (D(xA)\5(x,\)) 



+ 



g.2-p*,Q/3 



dx a d X0 G*'°(x,x) 



ci.«(r(x,4)\r(x,J)). 



From (5.2) and the boundedness of Gtl ' 7 , one immediately obtains 



(5.5) 

Now, for < < 1, 
(5.6) HZ,*'^" 



l z w 6 'lr»(fl(«,4)\E(*,i)) - ( £3 



62 \\C^(T(x,4)\T(x^)) 



d Xa d Xj3 G*>°(x,x) 



Ci.-'(r(x,4)\r(x,|)) 



O (e^) . 



Expanding F e in powers of e, one notices that the terms of order — 1 and in e cancel and 
that for all x £ D(x, 4)\D (x, |) 

F%x) =s\v s - A*(y)X7 x ( X ^ a (mx a G*'°(x,x)) 



+ V x -A*(y)V y [T*^{y)d~ Xa d~ XB G* 



[x, x 



Vy-A*(y)V £ (T*> a P(y)d Xa d Sg G*>°( 



x, x 



X, 



V* • A*(y)Vx r*^(y)^ Q d-G*' (x,x 



This demonstrates that F £ is (e) in L d+S (D (x, 4) \ D (x, \)). 

The source term Fi, due to the boundary layer deserves more attention. The expansion 
of F^ in powers of e is quite heavy. For the first term in the r.h.s. of (5.3) we get 



eV x -A*(-) V, Gw' 1 ' 7 ( - d*G*>°(x, x) 



(5.7a) 
(5.7b) 
(5.7c) 
(5.7d) 



W ( I ) GJ 1 ' 7 



+ A 



u yp^bi 



d Xa d x G*>°(x,x) 



+ A*<* ( \ ) S^G*' 1 ' 7 ( ? ) flk,^ G*>°(x,x) 



d Xa dx l3 dx- l G*'°(x,x) 



We write the second term in the r.h.s. of (5.3) as a sum of three terms 



e ^V 5 • A* - Vx X 



*,7 



C 



dx J dx v G* , ° (x, x 



T° (x, x) + eT 1 (x, x) + e 2 T 2 (x, x) , 
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with at order e 1 



(5.8a) T°(x,x) :-- 

(5.8b) 

(5.8c) 

(5.8d) 



+ 



dx-,dx v G*'°(x, x) 



dy a (A*^(y)dy^(y)) G* b ?«(y)] d^G*'\x, x) 

d ya (A*^(y) X ^(y)) dy G*f« (y) 
+ \A*^(y) X ^(y)d ya d y ^(y) 



+ 



dx-ydxr,G*'°(x,x) 



dx,d Xv G*'°(x,x), 



at order e 1 



(5.8e) T 1 (x,x) :-- 

(5.8f) 

(5.8g) 

(5.8h) 

and at order e 2 



A*^(y)dy 0X ^(y)Gl^(y) 
+ 



d Xa d^d $TI G*V(x.j) 



A*^(y) X ^(y)dy Gl^(y)\ I J J <),J),, <),..( r»li, r) 

d ya (A*^(y) X ^{y)) G^\y)] (f\ 0,. i),.i),,(r »{.r.x) 



+ 



+ 



A*^(y)x*' 7 (y)dy a G* b ^(y)] r)d^d~ Xj d~ Xn G*'\x,x) 



(5.8i) T 2 (x,x) := [A*^{y)x^{y)G* b ^{y)] f^j d Xa d Xfi dx^dx v G*' (x,x) 



We intend to show that all these terms are O (e K ), with k > 0, in L d+S (D (x, 4) \ D (x, \)). 
This seems tricky for some terms in the expansion above. Indeed, (5.7a) and ( 5.8a[ ) are of 



order 0(1), as we do not know more than G^l £ L°°(£l n ). However, the sum of (5.7a) 



and (5.8a) cancels: 



d ia (A*^(y)dy^(y)) G* b ^(y)} ( ^ W^G*' (*,*) 



d ya (A^{y)dy vX ^{y)) G* b ^(y)] (f) d Xfj d^G*>°(x,x) 



" )<% ltV ~ G*>\x,x). 



For the remaining terms, we use either (2.12 ), if at least one derivative of G* b } ,r! is involved, 



or (3.12 



if not. Therefore we proceed in the same manner for (5.7b ), (|5.7c ), (5.8b), (5.8c) 



(5.8d), (5.8f), (5.8h) on the one hand, and (5.7d), (5.8e), (5.8g), (5.8i) on the other hand 
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Let us estimate ( |57Fb| > in L d+S (D (x, 4) \ D (x, {)): by ( jOg) with k' = 1, 



< 



- , yp u \ s i ' ' ' 



^y G bi 



L^(D{xA)\D{x,\)) 



D{x,A)\D{x,\) 
8 

< c ' 



dx 



sup 



(z 1: ...z d _ 



dt 



< C 



sup 



(zi,...z d _ 



V^'^lMlzi,...^!, 



2 \ 2 

dt 



and 



< <?£2 



sup 



V^ 1 " (M (zi, . . . z d _i,t)) 



dx a d Xj G* ,0 (x, x 



L°°(D(x,4)\D(x,±)) 



0(1) 



from which we get, by interpolation, 



L d + 6 (D{x,4)\D(x,\)) 



O £2(^+*) 



We have gained a positive power of e for all corresponding terms listed above. In the other 
cases, the use of (3.12) deteriorates the exponent of e in the estimate, though it remains 
positive. Let us bound ( |5.7d| in L d+S [D[x,A) \D(x,\)): 

eA* 



L d + s (D(x,A)\D(x,l)) 



< Ce 

I D(x,4)\D(x,i) 
Hence we have shown 



G 



bl 



d+S \ d+<5 

dx I < 



<Z+<5 \ d+<5 

dtl <Cei 



o V e 



•MllL^(D( :E ,4)\D(x,i)) 
This bound, with (5.5), (5.6), the estimate on F £ and ( |5.4j ), boils down to 

( 5 - 9 ) ll V ^W S ' a 'IL<»(£>(! e ,3)\D( a i 1 |)) - Ce * 

with a positive, albeit small, k := min ^1 — v, 2 (d+5) ) • One can arwa y s take v sufficiently 
small so that 1 — v > 2 (d+g) • ^ s ^^bf' X ^ s U P t° ^ ne boundary <9f2 n , it follows from 
(3.17) and (5.9) that one can expand P e (x, x). For all x 6 -D (x, 3) \ Z? (x, |) , 



(5.10) 



P £ (x,x)-P°(x,x,-) -ePMx,x, X J -e 2 P 2 (x,x, X 



< Cs K , 



with at order e° 

[P°(x,x,y)] T := [A*^(y) + A*^(y) (y) + d^G*^ (y))] d if3 G*>°(x,x)n 
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at order e 1 

' A *^(y) ( X *.7 ( ^ ) + G *.i,7 (5) - 
+ A*^(y)d y X^(y) + A*^(y)dy, (x*'^y)Gl^(y)) (y)] d- Xp d^{x, x)n a , 



[P\x,x,y)] T 



and at order e 2 



[P 2 (x,i,y)]' :=A*>^(y) T*^(y)+ X ^(y)Gl^(y) d^d £v G*>°(x,x)n a . 



For y G Q n , y G <9f2 ra and e : = 



ey and x := are such that \x — x\ = 1. 



Applying now (5.10) with x and x defined like this, and rescaling the estimate in the 
variables y, y using the scaling properties of P £ and G*'° (see lemma 3.7), we finally have: 



Theorem 5.3. For all < k < there exists C K > 0, such that for all y £ £~l n ,o an d 
y G dSlnfl, 

(5.11) \P T (y, y) - A*(y)V y G*>°(y, y) ■ n - A*(y)V y ( x *(y) • V y G*>°(y, y)) ■ n 



A*my Kay) ■ V y G*'°(y,y)) ■ n - A*(y)V y (F*(y) • Vp*>°(y,y)) 



n 



-A*(y)V y (vTfflCtfm • V?G*'°(y, y)) ■ n 



< 



G K 



\y-y 



d-l+K ' 



6. Convergence towards a boundary layer tail 



6.1. The convergence proof. It follows from theorem 3.5 that the variational solution vu 



of (1.1 ) in fl n fl (see section 2.2 ) can be expressed by the mean of Poisson's kernel associated 



to —V • A(y)V- and O n . Hence, by theorem 5.3 for all y G f2 n , for all i G {1, . . . -/V}, 



(6.1) 



vbi,i{y) 



Pij{y,y)vo,j{y)dy = / Pji(y,y)v ,j(y)dy 

A* k (y)V y \G*>°(y, y) + X *(y) ■ V y G*'°(y, y) 



+ G^(y) ■ V y G*>»(y, y) + T*(y) ■ Vfi*'\y, y) 



Y*,0 / ~ 



+ X *(y)G* b f(y) ■ V 2 y G*'°(y,y) , . • n v^(y)dy + 



ki 



Ri(y,y)dy 



where for all y G dCl n , \Ri(y, y)\ < 



8Q n 

< C 



Ri(y,y)dy 



< C 



c 

\y-y\d-l + n 
1 



. The bound on the remainder Ri yields 



an„ \y - y\ d ~ 1+K 



dy 



-!X{0} \ Z ~ X ~ 



-dz < 



C 



(y ■ n) K 



1 (1 + 1* 



z/\2 



d-l+i 



-dz' < 



C 



y-n— too 



(y ■ n) K 



0. 



where z := M T y. It remains to handle the other terms. The boundary function vo is 
quasiperiodic along the boundary dQ n albeit not periodic. This suggests to use the fol- 
lowing lemma to take advantage of the ergodic properties related to the quasiperiodic 
setting. 



Lemma 6.1 ( |Sub74| theorem S.3). Let f = f(y) G M. be a quasiperiodic function 
Then, there exists M{f] G K such that for all tp G L l (R d ) , 

<p(y)f (-) dy M{f} [ ^{y)dy. 

./rod 



on 
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Let i, j G {1, . . . N} be fixed. We focus on the convergence when y ■ n — > oo of 



A *jk(y)^y G *U % V) ■ n vo,j(y)dy = I A*'^ p {y)dy p G*^{y, y)n a v 0J (y)dy 



3Q r , 



For all R > 0, for all y = y' + (y ■ n)n G fi n with y' G <9f2 n D 5(0, i?), taking z := M T y and 



^4*f p (y)dy fi G* k f (y, y)n a v ,j (y)dy 



A*^(M2)d^Glf{Mz, M(z', z d ))n a v 0J (Mz)dz 

M^-ixjO} 



^^(Mf -.0 



e e 



1 /M(^0) M^l)\ / , 



d li/3 G*f (M(z',0),M(0, 1)) n Q A*>^ ( M ( j,0 ) ) v 0J (M ( ^,0) ) dz' 



(6.2) 



+ 



d h/3 Glf{M(z',0),M(ez',l)) 



-d 1>p G*» {M(2',0),M(0 A))}n a A*f (M ( -,0 ) ) v 0lj (M ( ^,0) ) dz' . 



Let us show that the second term in (6.2) tends to when e — > 0, uniformly in z' G 
5(0, R) C For e sufficiently small such that eR < 1, 

/ [e? 1)/3 G*f (M(i', 0), M(ez', 1)) - (M(z', 0), M(0, 1)) 

n a A*f (M(-,o])v jM(-,0))dz 



< 



cf sup | ViV 2 G*'° (M(z', 0), M(u', 1)) | e\z'\dz' 



< CeR 

< Ce. 



sup 



C 



-dz 



-i w'es(o,i) (l + |F - u '|2) 



From the bound (3.5c), we get for all y G c?ri r 



|ViG*'°(y,y)| <C 



y ■ n 

\y - y\ c 



which shows that 



d hP G* k f (M(*',0),M(0,1)) 



< 



C 



(1 + \z< 



d ■ 

2U 



Hence d lt pG* k f (M(z' , 0), M(0, 1)) G L|, (M^ -1 ) and one can apply lemma 6.1 to get the 
convergence of the first term in (6.2): 



di,p G * k f (M(^,0),M(0, 1)) n a A*f ( M ( -,0 ) ) v ,AM ( -,0) ) dz 1 



£ 



e-s>0 



M 



dn n — > 
v •— ► A *jk^(.y) v o,j(y) n a J Jn n 
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d yp G kf(y^ n ) d y- 



The reasoning is identical for 



(6.3a) 
(6.3b) 



A *jk(y)^y X*kl{y)^yG*u{y,y) -nvQj(y)dy 



an, 



A *jk{y)^y G tiM(y)^y G *if(y^y) -nvoj(y)dy, 



as both terms involve just one derivative of G*'°. For (6.3b), we notice that y i— > £%„Gr^' ' 7 \y) 
is quasiperiodic on d$l n : we know from section 2.2 that there is a unique smooth V 7 = 
Vi\6, t) G R N , defined for 6 G T d and t > 0, such that for all y G <90 n , 

G* b i ln (y) = Gtf>\Mz) = v^m',o). 

The other terms in (6.1) involving strictly more than one derivative of G*'° tend to 
when e — > 0. Indeed, for all R > 0, for all y = y' + (y • n)n G Q n with y' G d£l n Pi i?(0, i?), 
again z := M T y and e := — - > 0, 

^jk^i^xl^i^dl^G^iy^^voji^dy 



dftn 



A*f(M~z) X l?(Mz)dl^G*/(Mz, M(z' , z d ))n a v^{M~z)d~z 



i*,0 



R^-ixjO} 



x A*f (M o)) xll 7 (M (~ o)) (M(^0),M(a',l)) 



dz 7 , 



which is easily shown to be of order 0(e). We proceed following the same method for 
the remaing terms. This demonstrates the convergence of the boundary layer towards a 
constant vector field v^, the boundary layer tail: for all y = y' + (y ■ n)n G Vt n with 
y' G d£l n n -B(0, R), 



(6.4) v bl (y) — »• w w 



=-/ 



dy a G°{n,y)dy 



M{A? a (y)v (y)np} 



+M 



{dy,(x*n T (y)A^(y)vo(y)n 1 } + M{d y ^ 



locally uniformly in y'. Moreover, (6.4) yields an explicit expression for in terms of the 
means Ai {•} on dVt n . 



1.2 



The 



6.2. The boundary layer tail does not depend on a. In the preceding section we have 
shown the convergence of defined in f2„ )a towards vj 00 when y ■ n — > oo. It remains to 
show that v^l°° is independent from a in order to complete the proof of theorem 
fact that n ^ M.Q d is crucial. To do so, we generalize lemma 6 in |GVMll] : 

f. Then, 

,a,oo ttoAT 



Proposition 6.2. Assume that n ^ K 

a G 

is Lipschitz continuous. 



v. 



bl 



G 



The proof in [GVMllJ relies on the small divisors assumption (1.12) and energy esti- 
mates. We, instead, have recourse to Poisson's integral formula and estimates on Poisson's 
kernel. Let a, a' G R and v := \a' — a\. We call G a (resp. P a ) the Green (resp. Poisson) 
kernel associated to — V- A (• + an) V- and Q, n = Q n fi. We define analogously G a ' and P a ' . 
We also have the aversions G*' a and G*' a ' corresponding to the transposed operator (see 
section 3.1). The following lemma is an adaptation of the results due to Avellaneda and 
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Lin (see |AL87] ). Kenig and Shen (see |KS11| section 2) and Gerard- Varet and Masmoudi 
(see |GVM| appendix A). 



Lemma 6.3. Let < /x < 1. There exists C > such that for all y, y G VL n , y ^ y, 



(6.5a) 
(6.5b) 



G a {y,y)-G a (y,y) 
G a (y,y)-G a '(y,y) 



< Cv 



~,d-2 1 

y-y\ 



ifd>3, 



(y ■ nf (y ■ nf 



< Cv 



\y-y 



d-2+2[i 



for all d>2. 



Proof. The ideas for the proofs are in large part taken from the reference above. For (6.5a 
we rely on a representation formula of G a (y, y) — G a (y, y): 



(6.6) G a (y,y)-G a (y,y) 



d 2 , a G a (y,y) (A^(y + na')-A^(y + na)) d^G a \y,y)dy. 



This formula, which is proven in |HK07| (see corollary 3.5) for the domain M. d , d > 3, is a 
consequence of Green's representation formula. The proof of Hofmann and Kim extends 



to the domain Q n . From (6.6) and (3.5d), one deduces that for all d > 3, y, y G Q n , y ^ y, 



G a (y,y)-G a (y,y) 



< Cv 



< Cv 



< Cv 



1 



1 



n n \y -y\ d 1 \y-y 



^dy 



1 1 

i = ^iT^id^dy 

\y ~ y ~ y\ \y\ 



l 



y-y\d 2 



\y-y\ y 



1 

d-1 ly^-l^ 



< Cv 



y-y\d 2- 



For (6.5b) we have recourse to the local boundary estimate (3.8). Assume that d > 3 



For given y, y G Q n , y ^ y, we first establish the bound 



(6.7) 



G a (y,y)-G a (y,y) 



< Cv 



(y ■ n) 



Let r := \y — y\. We distinguish between two cases. If y • n > |, then (6.7) follows directly 
from (6.5a). Assume that y ■ n < § and let y G dtt n such that y ■ n = \y — y\. Then, from 



(3.1) it comes that G a (-,y) — G a (-,y) satisfies 
(6.8) 

' -V y ■ A(- + na)V y (c a (•, y) - G a ' (•, y)) = V • [A{- + na) - A{- + no!)] VG a '(;y) 

in D{y,t) 

G a (-,y)-G a ' (-,y) =0 in T (y, |) 
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Applying a rescaled version of (3.8), one gets using (6.5a) and (3.5d), 

G a (y, y) - G a ' (y, y)\ = \ G a (y, y) - G a ' (y, y) - (G a ' (y, y) - G a (y, y)) 

\y- n \^n 



< 



G a (-,y)-G a (;y) 



< c 



d 

X 2 



C^(D(y,r.)) T» 

G a (;y)-G a '(-,y) 



+■ 



,d+5 



[A(- + na)-A(- + na')] VG a '(-,y) 



L*f*(D( S ,§)) 



< c 



< Cv 



G a (;y)-G a '(;y) 



+ rv 



VG ffl (-,y) 



y*D(y,i)\y-y\ d 2 ' ' y^D(i,i)\y-y\ 
\y ■ 

r d—2+(j, ' 



< Cv 



i'»(B(v,i)) 



\y ■ 

\y ■ n>\^ 



as for all y G D (y, |), |y — y\ > | , This shows (6.7). We now turn to the proof of 
(6.5b) itself. If y ■ n > |, then (6.5b) follows directly from (6.7). Assume that y ■ n < | 
and let y G <9f2 n such that y • n = |y — y\. Applying a rescaled version of (3.8) to 
G*' a (-,y)-G*' a ' (-,y) satisfying Q with A (resp. G a ', y, y) replaced by A* (resp. G*' a ', 
y, y), one gets using ( |3.5e[ ), and for all y G D (y, §) |y - y| > §, 



G a (y,y)-G a '(y,y) = G*' a (y, y) - G*' a '(y, y) 



< C 

< Cv 

< Cv 

< Cv 

< Cv 



G*> a (;y)-G*< a '(;y) 



L°°(D(y^)) 



+ rv 



VG*' Q (-,y) 



sup 



sup 



G a (y,y)-G a '(y,y) 



+ r sup 



y • n 



|y • 



(y ■ n] 



1-M 



(y • ra)^ (y • n) 

y^Z^ly-yl'-^ly-yl 1 -^ 



\y ■ n\ 



(y • 

(y • n) M (y • n) M 



Sl J P l 7 )_ U |d-2 + A1 +r 



|y • 



on condition that y ■ n < \y — y\ for all y € -D (y, 0. If the latter is not true, i.e. if there 
is y G D (y, |) such that y ■ n > \y — y\ > |, then we apply the same reasoning as above 
with A (resp. G a ', y, y) replaced by A* (resp. G* ,a ', y, y) to get 



(6.9) 



G a (y,y)-G a '(y,y) = G*> a (y,y) - G*> a ' (y,y) 



< Cv 



(y ■ nY 
\y-y\ d ~ 2+ ^ 



and deduce (6.5b) from (6.9) and y • n >£. The two-dimensional bound follows from the 
three-dimensional one as explained in |GVM| and [A L87| . □ 
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for £ G Z d , i?£ solving 



Let us notice that proposition 6.2 implies that for all a, a' G M, = Indeed, 



-V • A(y)Vvl = 0, y-n-a + £- n>0 



v 



bl 



satisfies, by periodicity of the coefficients and of the Dirichlet data, «£(•) := v b i{- + £)• 
vlf°. As n £ RQ d , the set {£ • n, £ G Z d } is dense in R. The 



Hence, f; 



i bl 



a— £-n,oc 



independence of v , ! from a follows now from the continuity of a \— > v. 



The rest of this section is devoted to the proof of proposition 6.2 Let as before a, a' G 



and v := \a' — a\. Let vu be the solution of (1.1) and v' b , the solution of (1.1) in the domain 



£l n , a ' instead of $l n , a - Then vu = v^(- — an) (resp. v^i = vf L (- — a'n)) where vfa (resp. v^) 
solves 



(6.10) 



-V • A(y + an)Vv^ = 0, 



y ■ n > 



v (y + an), yn = 



(resp. (6.10) with a replaced by a'). Take ip G (R), compactly supported in [—1,1] 
such that < (p < 1 and ip = 1 on [—5, |] • Then := — </j(y • n)vo(y + no) (resp. 
:= — <p{y • n)vo{y + na')) solves 



(6.11) 



-V • A(y + an)Vv bl = V • A(y + na)V (</j(y • n)vo(y + na)) , y • n > 



0. 



y ■ n = 



(resp. (6.11 ) with a replaced by a'). One important point is that the source term in (6.11 ) 
is compactly supported in the direction normal to the boundary. We now estimate 

<i(y) - vubj) = Vbiiv) - "hi (y) + <p(v ■ n ) [Mv + na ) - My + na')] 

for all t/ £ O n . We have, 

v a bl (y)-vu(y) 

G a (y, y)V • A(y + na)V (<p(y • n)t; (y + na)) ay 



(6.12a) 
(6.12b) 
(6.12c) 



G a '(y, y)V • A(y + na')V (p(y • n)v (y + na')) dy 
G a (y,y) - G a ' (y,y)jv ■ A(y + na)V {(p(y ■ n)v (y + na)) dy 
+ / G a '(y, y)V • [A(y + na) - A(y + na')] V (p(y • n)v (y + na)) dy 
+ [ G a '(y,y)V-A(y + na')V{<p(y ■ n ) \yo{y + na) - v (y + na')] ) dy. 
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We analyse the terms in r.h.s. separately. The first, (6.12a) deserves more attention. By 
estimate (6.5b) and the usual change of variables z = M 1 y, for y ■ n > 1, 



G a (y, y) - G a (y, y)J V • A{y + no) V (cp(y ■ n)v (y + na)) dy 
(y ■ nY (y • nY 



< Cv 



n n \y-y\ d ~ 2 + 2 » 



< Cv- 



| V • A(y + na)V (<p(y ■ n)v {y + na))\ dy 



dz' 



(y ■ n — l) 2 + \z'\ 2 



d-2+2^ 
2 



(y ■ n 



— l + U I d~2 + 2n 



du'. 



We need, 2u > 1 for the integral to be convergent, and u < 1 for the r.h.s. to be bounded 



when y ■ n — > oo. We now work with such a /j,. For (6.12b), we rely on (3.5c): for y ■ n > 1 



G a ' (y, y)V ■ [A{y + na) - A(y + na')] V (<p(y ■ n)v (y + na)) dy 

< Cv 



1 



- 1 fl + |n'| 2 l2 



-du' 



which is a convergent integral. For (6.12c), we argue analogously. We end with 



vti(y) - «« (y) < ««(*/) - «S' (y) + <p(y -n)Mv + na ) -Mv + na ')\ < 

for y ■ n > 1, which proves proposition |6.2| letting y • n — > oo. This concludes the proof of 
theorem 11.21 

7. Almost arbitrarily slow convergence 

We exhibit examples in dimension d = 2 showing that in general the convergence of vu 
towards its boundary layer tail vff can be nearly arbitrarily slow. Let us, for the rest of 
this section, focus on the case when n ^ MQ 2 . We take d = 2, N = 1, A = I 2 and study 
the unique variational solution v of 

-A z v =0, z 2 > 

v(z) = «d(Nzi), 22 = ' 

where as usual N £ E 2 is the first column vector of an orthogonal matrix M sending e 2 



(7.1) 



N • Vg 

d t 



on n. From theorem 2.2 we know that v = v(zi,z 2 ) = V(Nzx,Z2), with V = V(9,t) G 
(9,t) G T 2 x M + , solving 

2 

V= 0, i>0 

^(M)= «a(0), * = o 

Expanding i n Fourier series yields for all 9 G T 2 , 

Uo (0) = ^ tJoCOe*'**, 

where G l 2 (Z;R). From ( |2~6] ) , in particular d t V G L 2 (T 2 

(0, t) G T 2 x 



it comes for all 



V(9,t) = Y,vo(Oe- 27Tmit e 2i ^' 
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Parceval's equality 

(7-2) \\V(9,t)-e (0)\\ 2 = £ \vo(0\ 2 e- 4lSU 



|2 

€6Z2\{0} 



together with Lebesgue's dominated convergence theorem prove that 

\\V(9,t)- vo^Wl^^O. 

As vq is a C°° function, its Fourier coefficients (vo(^))^ go to zero when |£| — > oo, faster 
than any negative power of |£|. It follows from this, for all a G N 2 , 

(7.3) ||^ (mt )_,r (o))|| 2 2(T2) ^0. 

Using Sobolev's injections, one notices that ( |7.3[ ) proves again the convergence of v towards 
the boundary layer tail := vq(0). 

Assume for a moment that n satisfies the small divisors condition (1.12). Let us come 
back to (F7721). For all m G N, 



\V(e,t)-v (0)\\ 2 L2{T2) =t- m Yl \vo(0\ 2 t m e- 47Tmit 



?ez2\{o} 

|N-e| 



\ ■ < i'» 1 

5ez2\{0} 

<Ct' m \vo(0\ 2 \^ 2+T)m {\^-m m e- 4nmit 
«ez2\{o} 

(7.4) <cr m \MO\ 2 \Z\ {2+T)m <c m t- m , 

H& 2 \{o} 

the function t H> (|N • £\t) m e - 47r l N 'SI* being bounded on E + . We have shown that V(-,t) 
converges to vq(0) in I? (T 2 ), faster than every negative power of t. 

Assume now that n ^ IRQ 2 does not verify (1.12) and let I > 0. Hence there are points 
of the lattice Z 2 , except 0, which are as close to the line N • y = as we wish. The sum 
in the r.h.s. of (7.2) does not necessarily keep the trace of the exponential behaviour of 
its terms. We show that the convergence is at least as slow as the convergence of 1 1— > t~ l 
towards at 00. In fact we aim at proving: 



Theorem 7.1. Assume that n ^ IRQ does not satisfy (1.12). 

Then, there exists a smooth vq and a strictly increasing sequence (£m)m>i of positive real 
numbers, tending to 00, such that for all a G N 2 7 for all M £ N \ {0} , 

\\d$(v(e,t M )-v (o))\\ L2{T2) >tj*, 

where V is the solution of ( |7.1| ) associated to vq. 

Let us insist on the fact that theorem 7.1 holds for any n ^ IRQ 2 , which does not satisfy 
the small divisors assumption. The idea of the proof is to choose a family {yo{^j) ^ whose 

support, that is the set of subscripts ( £ Z 2 such that Uq(£) ^ 0, is sufficiently close to 
the line N ■ y = 0. We now construct a suitable sequence (£m)m>i of 1? \ {0}. The small 
divisors assumptions being not verified, for all M G N\ {0}, there exists £ G 1? \ {0} such 
that 

(7-5) |N-£l<^ir M 
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One can construct (£m) recursively: 

fi := argmin |£|, £ 2 := argmin |£|, £ M : = argmin 
£ez 2 \{o} ?ez 2 \{0} «£Z 2 \{0} 

IN-^^ICr 1 I6l>l?i|+1 \^m\>\^m-i\+1 



|n-«k|i?|- 2 |N-ei<i|«|- 



M 



where argmin stands for a minimizor. The existence of a minimizor is ensured by (7.5). 
For £1 the reasoning is straightforward. Let us sketch the proof of the existence of £2, 
which immediately applies, with minor modifications, for all £m- We assume that for all 
|e| > 161 + 1, |N-el > \\t\~ 2 - Then, for all M > 2, |N ■ £| > ^\~ 2 > ^|e|" M - According 
to (Tty , this shows that there exists £ G I? \ {0}, |£| < |£i| + 1, such that, for all M > 1, 
|N • £| < ^|£| -M . For triis £ / 0, N • £ = 0, which is incompatible with n ^ IRQ 2 . We have 
thus built a sequence (£m)m>i of vectors of Z 2 \ {0} satisfying: 

(1) (\Cm\) m is strictly increasing; 

(2) for all AT > 1, |£m| > M; 

(3) for all M > 1, |N • 6/| < jyl&lf l - ^- 

We now come to the construction of Vo keeping in mind that vq has to be smooth. For 
all £ £Z 2 , we define 

0, if £#6/, -6/ for allM> 1 



M-% M \- M \ if£ = 6i/or -e 



M 



Thanks to the construction of the sequence (£m)m, f° r ah m G N, wo(£) = O (|£| m ). Thus 
u defined like this is a C°° (T 2 ) function and V = V(9, t) defined by for all 9 G T 2 , for all 
t > 0, 

y(0,t) := 2 ^ M-'KMr^e" 2 ^-^ 1 ' cos(2tt£ m • 0) 

M>1 

is a smooth solution to ( |7.1[ ). 



For all M > 1, let ijy/ := I— ■ The final step in the proof of theorem 



7.1 



is to 

estimate \\dg (V(9,t]\j) — vo(0))\\ L2 ^ 2 ^ f° r a S N 2 and M > 1. Of course, in our example 
^o(0) = 0. Yet one is free to modify this coefficient without changing anything to the 
nature of the problem. Let a G N 2 . By Lebegue's dominated convergence theorem, for all 
(6>,t)GT 2 x 



+ j 



d%V(9,t) = 2^2 M-'|e M r M/ (2^)l Q l^ / e- 2 " |N '^ |i cos (|a|) (2< 

A/>1 

and Parceval's inequality yields 

(7.6) \\d%V(9,t)\\l 2{j2) = 2 £ ikT- 2 '|a/r 2A// (2vr) 2 l«l|^| 2 e- 4 -l N -^l*. 

Due to our choice of sequence (£m)m 5 m particular because of the second property of (£m) 
listed above, there exists M<M > 1, such that for all M > M* ', £m,i ^ and £ M ,2 ^ 0, 
where £ M = (£m,1,£m )2 ); for all M > M<®, \^\ = IQJIQ2I > L Therefore, (^ yields 
for all t G E+, 

||WM)||| 2(T2) >2 £ M-^lafl-^^vr) 2 !-!!^^- 4 ^-!* 

M>M(°) 

>2 Yl M- 2l \C M \~ 2Ml e- inmMlt 

M>M(°) 

>^r* Y ( 47Tt Ye-^>* 

^Am\zm\ m ) 



M>A/(°) 
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and for all M > M(°\ 
Theorem 



7.1 



\d%V{0,t M )\\ L 2 (T2) > V2(^) l tJ, which proves theorem [t7T| 



prevents V(9, t) from decaying fast towards when t — > oo: in particu- 
lar, (7.4) is impossible, in general, when n does not meet the small divisors assumption. 

Let us turn to estimates in 



However, theorem |1.3| cannot be deduced from theorem |7.1| 
L°° norm for v. Uniformity in the tangential variable 9 is replaced by local uniformity in 
z\. Let R > 0. We slightly modify (v (^))^. As for all \zi\ < R, 



27t|£a/-N||z 1 | < %t\Zm\ 



-M 



< 



2vri? 



< 



7T 



M ~ M " 4' 

for all M > sufficiently large, depending on R, we consider vq defined by 

0, if£^M, -Cm for allM>l 

0, if £ = £ M or - £ M for 1 < M < AfW . 

M-'ICmI-^, if e = £m or - 6/ for M > 
Then, for all |zi| < R, for all t G 



«o(0 ■-- 



v(zi,t) 



> 



2 £ M 

M>M(!) 

2 



-/ 



ICm| 



-MZ p -27r|N-C M |t 



COs(27T^Af • Nzi) 



(2vr) z 



E 

M>M(!) 



2vrt 



M\Zm 



M 



which demonstrates theorem 1.3 when evaluated in t = tM for M > 
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